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AVANT-PROPOS 



Bien que le developpement et la disponibilite de calculateurs aujourd'hui rendent 
indispensable qu'un ingenieur soit famitier avec les methodes numeriques de resolution de 
problemes de transfert de chaleur, !a comprehension des mecanismes entrainant ce 
transfer! restent la voie obligatoire a emprunter afin d'arriver a ce but. 

Le present ouvrage voudrait contribuer a repondre a la necessity de faire apprecier a 
I'etudiant les mecanismes qui commandent ce phenomene assez commun d’echange 
d’energie qu'est le transfert thermique. II est le fruit d'un support de cours enseigne a I'Ecole 
Nationale d'lngenieurs, a [’university d’Annaba et au Centre Universitaire de Guelma. II 
comporte principalement cinq chapitres dont les trois premiers traitent chacun d'un mode 
particulier de transfert de chaleur (conduction, convection et rayonnement) tandis que le 
quatrieme traite de I’echange de chaleur avec changement de phase. Les echangeurs de 
chaleur constituant le dernier chapitre de I'ouvrage, sont une application directe des 
differents phenomenes etudies qu'on rencontre dans de nombreux secteurs de I'industrie. 

L'auteur essaye aussi de repondre a une exigence generalement consideree comme 
une caracteristique essentielle dans un ouvrage destine a des etudiants en cycle 
universitaire: la disponibilite d'exercices resolus. Aussi, tous les chapitres ont ete illustres 
d'exemples traites dans le detail. En plus, des questions de comprehension ainsi que des 
problemes neprysentatifs de situations concretes ont yte inclus a la fin de chaque chapitre. 
Nyanmoins, le lecteur pourrait noter I'absence de programmes de calcul specifiques. En 
effet, sachant qu'i! existe plusieurs manieres d'approcher la solution d'un probleme de 
transfert thermique et assumant que ce dernier possede les principes de base lui permettant 
de rediger un programme de calcul quand I'algorithme est connu, l’auteur s'abstint d'en 
presenter et s'est tenu au traitement d’exemples representatifs des techniques de resolution 
les plus utilises. Le lecteur interesse pourra trouver un grand nombre de logiciels traitant 
divers problemes de transfer! de chaleur en consultant la litterature specialisee. 



L'ouvrage est enfin complete par une serie d'annexes rassemblant les principales 
donnees thermophysiques necessaires a la resolution d'un probleme thermique. II est cloture 
par des tables de conversion entre le systeme d'unites intemationales (SI) et les unites 
anglo-americaines encore tres utilisees. 
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NOTATIONS 



Sauf indication contraire, les notations utilisees sont : 



1 - Caracteres latins 



a 

c 

c o 

C\,C 2 

4 

c v 

d,D 

e 

E 

4 

G 

Ec 

h 

H 

1 

J 

k 

L 

4 

4 

m 

M 

M, 

71 

P 

P 

<Iv 

r 

R 

s 

5 

i 



Diffusivite thermique 

Chaleur massique ou vitesse de propagation des ondes electromagnetiques 
Vitesse de propagation des ondes electromagnetiques dans le vide 
Constantes de la loi de Planck 

Chaleur specifique du liquide a la temperature de saturation 

Chaieur massique a pression constante 

Constante 

Diametres 

Epajsseur 

Eclairement 

Facteur de forme 

Vitesse de masse 
Constante 

Coefficient d'echange de chaleur par convection ou constante de Planck 
Enthalpie 

Intensite d'un rayonnement 
Radiosite 

Coefficient d'echange de chaleur par conduction 

Longueur ou luminance d'un rayonnement 

Chaleur latente de condensation 

Chaleur latente de vaporisation 

Debit massique 

Emittance 

Emittance monochromatique 

Indice de refraction ou normale a la surface 

Pression 

Pression ou peri metre 
Quantite de chaleur 
Debit volumique 
Coordonnee ou rayon 
Constante des gaz parfaits 
Constante 
Surface 
Temps 





Notations 



T : Temperature 

T ce : Temperature d'entree du fluide chaud 

T cs : Temperature de sortie du fluide chaud 

Tf e : Temperature d'entree du fluide froid 

Tf S : Temperature de sortie du fluide froid 

: Temperature de la phase gazeuse 
T p : Temperature a la paroi 

T sal : Temperature de saturation 

: Temperature loin de la paroi 
U : Vitesse 

u : Composante de vitesse/x 

v : Composante de vitesse/y 

V : Volume 

W : Travail 

w : Composante de vitesse/z 
x,y,z : Coordonn6es 



2- Caract6res grecs 

a : Absorptivite 

a k : Absorptivite monochromatique 

p : Coefficient de dilatation volumique 

5 : Epaisseur ou epaisseur de la couche limite 

s : Emissivite ou efficacite d’une ailette 

X : Longueur d'onde 

ji : Viscosite dynamique 

v : Viscosite cinematique ou frequence de radiation 

p : Masse volumique ou reflectivite totale 

p k : Reflectivite monochromatique 

a : Constante de Stefan-Boltzman ou tension de surface 

t : Contrainte de cisaillement dans un fluide ou transmittivite totale 

<p : Densite de flux thermique 
: Flux thermique 
¥ : Angle 

ft : Angle solide 



3- Nombres sans dimensions 

Bi : Nombre de Biot 

Gr : Nombre de Grashoff 

Nu : Nombre de Nusselt 

Pr : Nombre de Prandtl 

Ret : Nombre de Rayleigh 

Re : Nombre de Reynolds 




SOMMAIRE 



Notations utilisges 
Sommaire 



CHAPITRE : 1 

INTRODUCTION AUX TRANSFERTS THERMIQUES 



1.1 


: Introduction 


1 


1.2 


: Rappels de thermodynamique 


1 




1.2.1 


Le premier principe de la thermodynamique 


1 




1.2.2 


Le second principe de la thermodynamique 


2 


1.3 


: Diff6 rents modes de transfert de chaleur 


2 




1.3.1 


La conduction 


2 




1.3.2 


La convection 


2 




1.3.3 


Le rayonnement 


2 


1.4 


: Conclusion 


3 






Chaoitre : 2 








TRANSFERT DE CHALEUR PAR CONDUCTION 




2.1 


: Concept de la conduction 


4 


2.2 


: Lois fondamentales de la conduction 


4 




2.2.1 


Definitions 


4 




2.2.2 


Loi £lementaire de Fourier 


5 




2.2.3 


Equation generate de la conduction 


7 


2.3 


: Conditions aux li mites 


10 


2.4 


: Regime permanent-conduction stationnaire 


10 




2.4.1 


Le mur 


10 


2.5 


: Convection de la chaleur des surfaces 


12 


2.6 


: Rayonnement de la chaleur des surfaces 


13 


2.7 


: Exemples 


13 




2.7.1 


Mur simple en contact avee deux fluides 


13 




2.7.2 


Mur compose 


14 




2.7.3 


Cylindre creux a surfaces laterales isothermes 


16 




2.7.4 


Sphere creuse a surfaces isothermes 


20 



Sommaire 



IV 



2.8 : Probleme de conduction unidimensionnelle 

avec dfeperditions h travers les surfaces laterales 23 

2.8.1 : Equation generale 23 

2.8.2 : Solution elementaire 

pour les ailettes unidimensionneiles 24 

2.8.3 : Efficacite d'une ailette 25 

2.9 : Conduction bi-dimensionnelle, stationnaire, sans 

sources - met h ode de resolution par separation des variables 28 

2.10 : Conduction unidimensionnelle en regime variable, 

sans sources: methode de la transformee de Laplace 32 

2.10.1 : Generates sur la transformee de Laplace 32 

2.10.2 : Application de la transformee de Laplace a (’equation 

unidimensionnelle de la chaleur en regime variable 33 

2.11 : Methodes numeriques de resolution de liquation de la chaleur 35 

2.11.1 : Generates 35 

2.11.2 : Etablissement des equations aux 

differences finies en regime permanent 35 

2.11.3 : Conditions aux limites 37 

2.11.4 Methode de resolution: 

Etude d’un cas pratique 38 

2.1 2 - Regime variable 44 

2.12.1 : Generates 44 

2.12.2 Methode explicite 45 

2.12.3 : Methode im pi icite 46 

2.12.4 : Methode generale 46 

Exercices 46 



Chapitre : 3 

TRANSFERT DE CHALEUR PAR CONVECTION 

3.1 Introduction 57 

3.2 Equations generates de I'hydrodynamique 57 

3.2.1 : Generates - Definitions 57 

3.2.2 : Equations de conservation 58 

3.2.3 : Le concept de la couche limite et ses hypotheses 61 

3.2.4 : Regime laminaire et regime turbulent 63 

3.2.5 : Nombres sans dimensions 66 

3.3 : Donnees empiriques de la convection 67 

3.3.1 : Utilisation du nombre de Nusselt 67 

3.3.2 : Determination du Nusselt 68 

3.4 : Methodes devaluation de h 69 



Sommaire 



3.5 : Analyse dimensionnelle 

3.5.1 : Determination du nombre de groupes adimensionnels: 

Application du theoreme de Buckingham ou du n 

3.5.2 : Exemple d’application du theoreme de Buckingham 

3.6 : Formules empiriques couramment utilisees 

3.6.1 : Convection forcee 

3.6.2 : Convection libre ou naturelle sur parois isothermes 
Exercices 



69 

70 
70 

72 

72 

77 

80 



Chapitre : 4 

TRANSFERT DE CHALEUR PAR RAYONNEMENT 



4.1 


: Introduction 


87 


4.2 


Rayonnements electromagn6tique et thermique 


88 




4.2.1 : Spectre des ondes electromagnetiques 


88 




4.2.2: Le rayonnement thermique 


68 


4.3 


: Classification des grandeurs physiques utilisees 


89 




4.3.1 : Classification suivant la 






composition spectrale du rayonnement 


89 




4.3.2 : Classification suivant la 






distribution spatiaie du rayonnement 


89 


4.4 


: Definition des grandeurs physiques utilisees 
4.4.1 : Grandeurs relatives a une 


89 




surface emettant un rayonnement 


89 




4.4.2 : Grandeurs relatives a une 






surface recevant un rayonnement 


92 


4.5 


: Lois du rayonnement thermique 


92 




4.5.1: Le corps noir 


92 




4.5.2 : La loi de Planck 


93 




4.5.3 : Lois de Wien 


94 




4.5.4 : Loi de Stefan-Boltzman 


95 




4.5.5 : Realisation pratique du corps noir 


95 


4.6 


: Emittance de corps reels 


96 


4.7 


: Notions de transmission, d’absorption et de reflexion 


98 


4.8 


: Le facte ur de forme et le flux net 

radiatif echang6 entre deux surfaces noires 


100 




4.8.1 : Le facteur de forme 

4.8.2 : Le flux net radiatif 


100 




echange par deux surfaces noires 


102 




4.8.3 : Representation des echanges radiatifs 






entre surfaces noires par analogie electrique 


103 




Sommaire 



VI 



4.9 : Le flux net radiatif echang£ par deux 

surfaces grises et sa representation 6lectrique 104 

4.9.1 : Le flux net radiatif 

echange entre deux surfaces grises 104 

4.9.2 : Representation analogique des 

echanges radiatifs entre surfaces grises 105 

Exercices 107 

Chapitre : 5 

TRANSFERT DE CHALEUR AVEC CHANGEMENT D’ETAT 

5.1 : Introduction 112 

5.2 : Transfert de chaleur par condensation 112 

5.2.1 : Introduction 112 

5.2.2 : Transfert de chaleur par condensation ie long de plaques 

et cylindres verticaux en regime laminaire 113 

5.2.3 : Transfert de chaleur par condensation le long de 

cylindres horizontaux en regime laminaire 116 

5.2.4 : T ransfert de chaleur par 

condensation en regime turbulent 119 

5.2.5 : Condensation dans les conduites horizontales 121 

5.3 : Transfert de chaleur avec Ebullition 122 

5.3.1 : Introduction 122 

5.3.2: Differents regimes d’ebullition 122 

5.3.3 : Caracteristiques a I’equilibre d'une bulle 123 

5.3.4 : Correlation pour une ebullition 

nuclEee en regime laminaire 124 

5.3.5 : Ebullition nucleee en convection forcee 125 

5.3.6: Flux maximal transfe re 128 

5.3.7: Ebullition par film 129 

5.3.8 : Cas particular de I’ebullition de I’eau 129 

Exercices 131 

Chapitre : 6 

ECHANGEURS DE CHALEUR 

6.1 : Introduction 134 

6.2 : Coefficient global d’Echange 134 

6.3 : Principaux types d’Echangeurs et leur classification 137 

6.3.1: Principaux types d’echangeurs 137 

6.3.2: Classification des echangeurs 138 





Sommaire 



vii 



6.4 : Distribution de temperature dans un dchangeur 

6.4.1 : Echangeurs a courants parallels 

6.4.2 : Echangeur a contre-courant 

6.5 : Evaluation des performances thermiques 

d’un echangeur en regime permanent 

6.5.1: But de revaluation 

6.5.2 : Hypotheses de calcul 

6.5.3 : Methode de la moyenne logarithmique, 

6.5.4 : Methode du nombre d’unites de transfert, NUT 

6.5.5 : Autres types d'echangeurs 

Exercices 



B1BLIOGRAPHIE 

ANNEXES 

A : Quelques dates 

B-1 : Proprietes de quelques solides a 20°C 

B-2 : Proprietes de I'eau sous pression de saturation 

B-3 : Proprietes de I'air £ la pression atmospherique 

C-1 : Emissivites totales normales de diverses surfaces 

C-2 : Emissivites totales hemispheriques 

D : Facteurs de conversion 



139 

139 

140 

140 

140 

141 
141 
145 
150 

150 



155 

157 

158 

159 

160 
161 
162 

164 

165 





Chapitre 1 : 

INTRODUCTION AUX TRANSFERTS THERMIQUES 



• Introduction 

• Rappels de thermodynamique 

• Differents modes de transfert de chaleur 

• Conclusion 



1.1 Introduction 

Le transfert de chaleur est I'un des modes les plus communs d’echange d'energie. II 
mtervient naturellement entre deux systemes des qu'existe entre eux une difference de 
temperature et cela quel que soit le milieu meme vide qui les separe. 

A la base de I’etude des transferts thermiques se trouvent les principes de quantity 
de chaleur et de difference de temperature definis par la thermodynamique dans ses 
principes memes. Toutefois, la thermodynamique classique traite de I’etat des systemes 
sous I’angle macroscopique et ne fait aucune hypothese concernant la structure de la 
matiere. Ainsi, elle s’attache aux etats d’equilibre et neglige les differents mecanismes qui y 
conduisent. L’etude de ceux-ci constitue la discipline qui nous interesse: les transferts 
thermiques. 



1.2 Rappels de thermodynamique 

1.2.1 Le premier principe de la thermodynamique 

Le premier principe est fondamentalement un principe de conservation. II regit les 
transformations mutuelles du travail et de la chaleur en etablissant un bilan. II peut 
s'enoncer comme suit: 

"Lorsqu'un systeme isole subit une serie de transformations qui le ramenent a un etat 
d'equilibre final (2) identique a I'etat d'equilibre initial (1) et au cours desquelles il n'echange 
avec I'exterieur que du travail et de la chaleur, le rapport du travail a cette quantite de 
chaleur est constant” 

\W + O] ‘ = Constante (T 1 ) 



W : Travail echangee entre le systeme et le milieu exterieur 
O : Chaleur echangee entre le systeme et le milieu exterieur 
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1.2.2 Le second principe de la thermodynamique 

Le premier principe de la thermodynamique pemnet de faire le bilan des echanges 
d'energie d'un systeme materiel subissant des transformations quelconques, mais il ne fait 
aucune difference entre les di verses formes sous lesquelles cette energie se trouve fournie 
au milieu exterieur. Or, il n'est pas indifferent de recueillir la diminution d'energie d'un 
systeme sous forme de travail ou de chaleur comme il n'est pas non plus sans importance 
que la quantite de chaleur eventuellement fournie le soit a une temperature faible ou a une 
temperature elevee. 

Le second principe permet de calculer la quantite de travail susceptible d'etre 
recueillie dans une transformation determinee. II permet aussi de preciser le sens dans 
lequel tendent reellement a evoluer les systemes (alors que le premier principe qui est 
essentiellement un principe d'equi valence, ne nous indique rien sur le s les systemes (alors 
que le premier principe qui est essentiellement un principe d'equivalence ne nous indique 
rlen sur le sens de la transformation). 



1 .3 Diff6 rents modes de transfer! de chaleur 

Le transfert de chaleur peut etre defini comme la transmission de I'energie d'une 
region a une autre sous I'influence d'une difference de temperature. 11 est regi par une 
combinaison de lois physiques. 

La literature traitant du transfert de chaleur recon na it essentiellement trois modes de 
transmission de la chaleur : la conduction, le convection et le rayonnement 



1.3.1 La conduction 

La conduction est definie comme etant le mode de transmission de la chaleur 
provoquSe par la difference de temperature entre deux regions d'un milieu solide, liquide ou 
gazeux ou encore entre deux milieux en contact physique. 



1.3.2 La convection 

Les phenomenes de convection interviennent dans la transmission de la chaleur 
chaque fois qu’un fluide se deplace par rapport d des elements fixes. Lorsque se produit au 
sein du fluide des courants dus simplement aux differences de densite resultant des 
gradients de temperature, on dit que la convection est naturelle ou libre. Par contre, si le 
mouvement du fluide est provoque par une pompe ou un ventilateur, le processus est appele 
convection forcee. 



1.3.3 Le rayonnement 

Le rayonnement est le mecanisme par lequel la chaleur se transmet d'un milieu a 
haute temperature vers un autre a basse temperature iorsque ces milieux sont separes dans 
I'espace. Ce mode de transfert ne necessite pas de support materiel et peut done s'effectuer 
dans le vide. En general, les sources de rayonnement sont des solides et le rayonnement se 
fait par la surface. 
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1.4 Conclusion 

La plupart des phenomenes etudies font apparartre I'intervention des trois modes de 
transmission de la chaleur en meme temps. Afin d'expliciter ceci, prenons le cas d'une 
chambre (of. figure 1.1) 

En effet, les trois modes de transfert de chaleur interviennent dans ce cas et sont 
represents par : 

■ La conduction a travers les murs, le parquet et le toit de la 
chambre, 

• La convection qui prend place grace au mouvement de fair a 
I'exterieur et a I'interieur de la chambre, 

• Le rayonnement de la lampe eclairant la chambre. 





Conduction 

Convection 






Rayonnement 



Figure 1.1: Illustration du phenomene de transfert de chaleur se produisant 
dans une chambre close 
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• Concept de la conduction 

• Loi fondamentale de la conduction 

• Regime permanent - conduction stationnaire 

• Convection de la chaleur des surfaces 

• Rayonnement de la chaleur des surfaces 

• Conduction unidimensionnelle avec deperdition aux surfaces laterales 

• Conduction bi-dimensionneile stationnaire - Separation des variables 

• Conduction unidimensionnelle en regime variable - Transformee de Laplace 

• Methode numerique de resolution de S’equation de la chaleur - methode des 
differences finies 



2.1 Concept de la conduction 

Le transfert de la chaleur par conduction est associe aux mouvements de 
translation, vibration et rotation des molecules et atomes constituant le milieu considere. 

Quand la temperature d'un corps croit, I'energie cinetique des molecules et atomes 
le constituant croit aussi entraTnant une augmentation du transfert de chaleur par 
conduction. 



2.2 Lois fondamentales de la conduction 
2.2.1 Definitions 

Considerons une plaque (D) de surface S. Soft dO la quantite de chaleur echangee 
entre la plaque et I'air ambiant pendant le temps dt . 




Figure 2.1; Plaque plane 
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On definit: 

• Le flux thermique comme la puissance echangee par la surface S de la plaque. II 
s exprime: 



dQ . , 

<D= — ; d>] = Watts 
dt 



( 2 . 1 ) 



• La densite de flux thermique comme la puissance echangee par une surface unite de la 
plaque. Elle s’exprime: 



<P 



dO 

Jdr 



; [rp] - Watts ! m 2 



( 2 . 2 ) 



• La source de chaleur par la puissance thermique qu’elle produit. Dans le cas d’une 
reaction chimique, elle s’exprime: 

O , = A 0 .e~ at (2.3) 

ou: A 0 et a sont des constantes. 



2.2.2 Loi elementaire de Fourier 

Les principes fondamentaux de la thermodynamique nous font savoir que: 

• i- energie est conservee en I’absence de source de chaleur, 

• La chaleur transmise passe toujours du corps chaud vers le corps froid. 

Si on considere la plaque (D): 




Figure 2.2: Plaque plane a travers laquelle se transmet la chaleur 



On peut verifier experimentalement que: 

0 = k~(T y -Z) (2.4) 

u 

avec: 

• ?i T 2 

• Q: la quantite de chaleur transferee a travers (D). 

• k : le facteur de proportionnalite appele conductivite thermique qui est une 

caracteristique du materiau. 
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Pour un element infinitesimal dS , la relation (2.4) s'ecrit: 



8T 

dQ = -kcIS- 

on 



(2.5) 



(.’equation (2.5) est attribute au mathematicien franpais Jean Baptiste Fourrier qui 
en 1822 enonga sa loi qui peut se traduire comme suit: ”En tout point d’un milieu isotrope, 
le densite de flux thermique instantanee est proportionnelle a la conductivity thermique k 
du milieu et au gradient de la temperature’’. 

Ceci nous conduit a la forme vectorielle de la loi de Fourier qui exprime la densite du 
flux thermique definie comme la quantity de chaleur transmise par unite de surface. Soit: 

( 3~-kGmdT ( 2 . 6 ) 



Remarque: 

Par convention p est compte positivement dans le sens d’ecoulement de la chaleur 
e’est-a-dire vers les temperatures decroissantes. GrcidT est un vecteur porte par le meme 
axe mais de sens contraire a <p, d’ou le signe negatif de la loi de Fourier. 



Exempts 2.1: 

Catculez le flux thermique ainsi que la densite du flux thermique a travers une plaque 
plane et homogdne dont Itepaisseur est de 50mm si elte est (a) en acier inoxydable 
(k - 1 6W fmK) de dimensions 3mx2m ou en (b) en bbton (k = 0 , 921 V ImK) de dimensions 
30mx2bn, Dans les deux cas, les temperatures aux surfaces de la plaque sont maintenues 
constantes et igales £ T pl = I00‘C et T pl = 90° C. 



Solution: 



(a) Plague en acier: 

• Q = ® -T pl )= 16x^(100 - 90) = 19200 Worn 




19 200 
3x2 



= 3200 W/m 2 



(b) Plague en byton: 

• Q = ® ffo, - 7 -J = °. 92 ^( 100 - = 1 10900 Wa tts 




1 10 400 
30x20 _ 



= 184 W / nr 



La densite du flux de chaleur montre la capacity d’un materiau a transmettre de la 
chaleur par conduction. De ce cas, on peut deduire que I’acier est nettement plus 
conducteur que le beton. 
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2.2.3 Equation generate de la conduction 

Dans un systeme de coordonnees cartesien (0,x,y,z), I’equation de la conduction 
de la chaleur ou simplement I’equation de la chaleur peut etre developpee en considerant 
un volume infinitesimal v et en effectuant le bilan thermique relatif a ce volume pendant le 
temps dt . 




Figure 2.3: Bilan thermique £ travers un volume infinitesimal retativement 
au systeme de coordonnees cartesien 

Relativement a la section perpendiculaire a I’axe des jc, les quantites de chaleur 
s’expriment grace a la loi de Fourier: 

d Q, = k (tfydz) 

dO x+dx = dO x + -~(dO x )cbc+... (serie de Taylor) 

( \ dT 3 \ i v dT\ 

42U* = ~k (dydz) — + — k [dyct) —\dx 

d\ dT 1 

j;[ T+ i^ dx r ,dz 

Le bilan thermique relativement a I'axe des x est: 

&T &T 8~T 

dQ x - dQ x+dz = ~k — ^ + k—dydz + k —rdxdydz 

d\ 3T 1 

dQ x ~dQ^=— k— dxdydz 



(2.7) 
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En exprimant les differentes quantites de chaleur mises en jeu relativement aux 
parois perpendiculaires aux axes y et z , et en developpant les equations de la meme 



<iQ,-dQ„,=- 

I 



3T 

* — 
oy 



dz 



on aura: 




dxdydz 


(2.8) 


dxdydz 


(2.9) 



Le bilan des quantites de chaleur transmises par conduction a travers le volume v 



est: 



d 



dx 



k £L ) + £*£l ] + £*£T 

dx) oy) dzj 



dxdydz 



Si par unite de temps et de volume, la quantity Q'{x,y,z,t) est fournie, alors la 
quantite de chaleur generee dans I’element sera: 

Q dxdydz (2.10) 

La chaleur transmise dans l’6lement de volume due a la conduction et la quantite 
generee ensemble servent a augmenter I'energie interne de (’element. Cet accroissement 
s’exprime: 

dT 

pc dxdydz-^ (2.11) 



Le bilan energ6tique final donne: 



dT 



dT) d 



P c = k TZ + T. 



r aT \ 

k 



dt 3t\ dx) dy \ dy) dz \ dz) 



d 
+ — 



df^ 



+ <? 



(2.12) 



qui est (’equation generate de la conduction en coordonnees cartesiennes. Elle est 
gdnerale et est valable pour les milieux h&erogenes anisotropes car: 



k = k(x,y,z,t) 

< P = p(x>y,z,t) 
c = c(x,y,zj ) 



Pour un corps isotrope et homogene, I’equation precedente devient: 

d t 

dt 

?f = a^T + V (2.14) 

dt pc 

ou: a = kl pc est le coefficient de diffusivite thermique. 



k fd z T d l T d 2 T 



— + ' 



pc\dx * dy 



i + ■ 



+ ■ 



Q_ 

pc 



(2.13) 
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2.2.3. 1 Equation generate de la conduction en coordonnGes 
cylindriques 




Sachant que: 

x = r cos(0) 

< y = r sin(0) 
z ~ z 

(’equation de la chaieur s’ecrit en coordonnees cylindriques: 

ST f S 2 T }_ST_ J_S^T_ S 2 'A Q_ 

St \ Sc 2 + r Sr + r 2 S9 l + Sz 1 J + pc 



(2.15) 



2.2.3.2 Equation ggngrale de la conduction en coordonnees 
sphgriques 




Sachant que: 

x = r cos (p) sin (^) 
^ y~r sin(p)sin W 
z-r cos(^) 



2 Transfert de chaleur par conduction 



JO 



(’equation de la chaleur s’ecrit en coordonnees spheriques: 



d 



£T_ 1 d'W 1 1 

dt a |r Br 2 + r 1 sin(^) By 



sin(^) 



BT 

By 



B 2 T\ O' 



+ ■ 



/" sin 



V) 



pc 



(2.16) 



2.3 Conditions aux limites 

L’equation generate de la conduction est une equation aux derivees partielles, 
lineaire de deuxieme ordre. Elle admet une infinite de solutions et ne peut avoir de sens 
physique que pour des conditions definies representant les causes qui detemninent 
revolution du phenomene; ce sont les conditions aux limites. 



2.4 Regime permanent - Conduction stationnaire 

En regime permanent (BTIBt-0) sans sources (Q r - 0), la repartition des 
temperatures dans un milieu isotrope et homogene obeit a (equation; 



B l T B 2 T 8 Z T 
Bx 1 + By 2 dz 1 



(2.17) 



Dans ce qui suit, nous ne nous interesserons pas a la resolution de [’equation dans 
le cas general mais considererons plutot des cas simples qui souvent sont les plus utilises 
dans la pratique. 



2.4.1 Le mur 

Le “mur” est un milieu conducteur homogene limite par deux plans paralleies infinis 
maintenus a une temperature uniforme. 




X 

Figure 2.6: Distribution de la temperature t travers un mur en regime permanent 



Les sources de chaleur assumees inexistantes en regime permanent, I'equation 
fondamentale de la chaleur notee (2.13) se reduit a: 

T d 2 T 
Bx 2 dx 2 



(2.18) 
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qui possede une solution de la forme: 

T(x) = C,x + C 2 

Les constantes Q et C 2 sont evaluees a partir des conditions aux limites: 



(2.19) 



T=T X 

T=Z 



pour 

pour 



x = 0 
x = L 



Appliquant ces conditions, on obtient: 

„ Z-T 
C, = T x et C, = 



dou: 



T - T 

T(x) = ±~±x + T, 



( 2 . 20 ) 



Le flux de chaleur est constant en regime permanent et la loi elementaire de Fourier 
permet de I’exprimer: 



<p = Q = ~kS 



_dT 

cbc 



( 2 . 21 ) 



ciT T 1 T 

— = 2 ... 1 est la pente de la droite representee en figure 2.6. 
dx L 



et done: 



<D = O = kS 



T x -T z 



Exemole 2.2: 

La density du flux thermique & travers un mur plan d'epaisseur 50mm est 70 W / nr . 
Calculez la difference de temperature aux surfaces du mur et les valeurs num&riques du 
gradient de temperature dans celui-ci si ce mur est en (a) en laiton (k = IOOW'7 mK), (b) 

granit (* = 2,5PF / mK) et (c) en bois (k = 0,23 W / mK) . 

Solution: 

(a) mur en laiton : 
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\GrcidT\ = 



^L\ 

dx' 



AT _ 0,035 
e ~ 0,05 



0,7 Kim 



(b) mur en qranit : 

• A7=^^ = l,4 £ 
2,5 

(c) mur en bois : 



Ar= 7^O05 = 1 

0,23 - — - 



• \GradT\ 



1,4 

0,05 



= 28 Kim 



» 



Gradl\ - 



15,2 

— — = 303,3 Kim 
0,05 — 



Le bois montre une bonne aptitude pour (Isolation suivi par le granit et enfin par le 
laiton. Ce dernier ne fait "descendre” la temperature que de 0,7 K par metre de materiau 
(compare aux 287" et 304, 3 A" du granit et du bois respective me nt). 



2.5 Convection de la chaleur des surfaces 




Figure 2.7: Distribution de la temperature a travers un film mince 
de fluide adherant a une paroi solide 

Bien que le probleme de la conduction ait pour objectifs principaux la determination 
de la distribution de la temperature dans les corps ainsi que celle du flux de chaleur, les 
conditions de surface sont importantes. 

Dans le cas ou un fluide s’ecoule sur une surface solide, la distribution de la 
temperature est celle illustree par la figure 2.7, Au sein du milieu fluide, le gradient de 
temperature est confine a une couche mince situee pres du mur alors que plus loin de cette 
paroi, la temperature peut etre consideree constante. 

Dans le but de simplifier la distribution de la temperature, son profil est re m place par 
une droite illustree en pointing sur la figure 2.7. Physiquement, ceci suppose qu'un film 
mince d’epaisseur & adhere au mur et c'est a travers ce dernier que (’evolution de 
temperature prend place, Le profil de cette derniere etant assume represente par une 
droite, on peut supposer que le transfert de chaleur se realise sous forme de conduction 
comme dans le cas d'un mur solide. Ceci simplifie enormement les mecanismes reels qui 
prennent place et permet d’ecrire: 
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0 = 



k 



Jluhie 



S' 




( 2 . 22 ) 



La quantite de chaieur O peut done etre determinee si S' est connue. Ce dernier 
parametre depend de plusieurs facteurs dont la vitesse du fluid©, 9a forme de la surface 
ainsi que sa structure. En pratique, on prefere travailler avec ie rapport {k nmd J S’) plutot 
que<5'. Si on note: 

h= k ft L ( 223 ) 

0 

h representant le coefficient de transfert de la chaleur par convection. 

On aboutit a une relation exprimant le flux de chaleur echange par convection entre 

un fluide et une paroi solide: 

A 

® = hS<l> = hS[T / -T m ) (2.24) 

qui est la loi de Newton. 



2.6 Rayonnement de la chaleur des surfaces 

Le flux de chaleur rayon ne par un milieu de surface S et de temperature T 
s’exprime grace a la loi de Stefa n-Boltzman: 

Q>=oSf* (2.25) 

ou: <j : Constante de Stefan-Boltzman. 



2.7 Exemples 

L’application des principes precedemment enonces peut etre mieux pergue par 
I’introduction d’exemples. 



2.7.1 Muren contact avec deux fluides : analogie 6lectrique 





Figure 2.7: Mur simple en contact avec deux fluides 
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Soit le mur d'epaisseur L et de conductibilite thermique k tel qu’illustre par la figure 
2.7. 



La conservation du flux thermique impose: 



qui peut s’ecrire: 



T \-T 2 


1 

i 


t ,- r 4 




1 


L 


i 




hy S 


kS 


hy S 




T \-T 2 




T 


T \~ T 4 


K. 


K* 


Kv2 


Kv i + Kd + R cvi 



(2.26) 



(2.27) 



ou par anaiogie avec I’electricite, R n] , R cd et R n . 2 sont appeles les resistances 
thermiques. 



La loi de Fourier ainsi que celie de Newton exprimant les quantites de chaleur 
transmises par conduction et par convection respectivement peuvent done £tre comparees 
a celie d’Ohm qui exprime la difference de potentiel entre deux points d’un reseau 
eiectrique en fonction de sa resistance et de I’intensite du courant qui le traverse. Ainsi: 



Ts 




ii 

i 

.JS 


Loi de Fourier 


i 

h S 




= Ty~T z : 


Loi de Newton 


1 


I 


Ii 




R 


I 


1 

to 

ll 


Loi d' Ohm 



2.7.2 Mur compost 



♦—2 > 


4 — * 


4 ► 


n 


T*' 


R 


K 


k 2 




R x 


r 2 


R, 









T,» 



Figure 2.8: Mur compost et son schema eietUK^e 
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La conservation du flux de chaleur implique: 

^( 7 ; - r 2 ) - Ar, - rj = Ar, |-(7; - r 4 ) 

t } -t 2 t 2 -t. t 3 -t a t x -t 2 t 2 -t 3 t.-t 4 

4 ~~ j4 _ ~ R l ~ R 2 - R, 

k x S k 2 S k 3 S 

^ _ Ty- T 2 + T 2-T i + T i -T A ^ T,-T 4 
R\ + Rr, + R 3 R 

avec: 

R = R] + R^ + R 3 
' 

i 

Le mur compose se comporte done comme un mur simple de resistance thermique 

R ~ 7?[ + Rq + /?3 . 




Soit le mur compost illustrd ci-dessus. En assumant une conduction 
unidimensionnelle et connaissant les temperatures aux parois T pl et T p2 ainsi que les 

conductibilites thermiques de ses diffe rentes sections, calculez le flux de chaleur par unitd 
de surface ct travers le mur. 

Donnees: 

|r vl = 1850*C ; kA^mW/wK ; kC = 901V/mK ; kE = 120W/mK 
[r,=25’C ; kB = 45W/mK ; kD = 60W/mK ; kF = \65W/mK 




( 2 . 28 ) 



( 2 . 29 ) 
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Solution: 

L'analogie existant entre ies grandeurs thermiques et electriques nous permet de 
representer le schema equivalent au mur considere. 



R: 



Ra 

T p \ * AA~ 



i-V\A 



Rc 



R D 

AA — i 



Rr 



Re 






pi 



Le flux s'exprime done: 



<!>=■ 



avec: 



T -T 

l p\ w 

R a + R^+R? 



0,1 



R ~ L * ~ 

J k A S A 180x0,01x0,5 



Rc 



Rc 



0,2 



Re = 



dou: 



k c S c 90x0,005x0,5 

L s 0,2 



= 0,111 K/W , A ? = 
= 0,888 A' /If , A r = 
= 0,666 A' 1 11’ , A, = 



k E S s 120x0,005x0,5 



(a h + R n ){R c + R e ) 

R — = 1 036A » 

*« A„+A 0 +A,- + A r 



r b 


0,2 


k ? S s 


_ 45x0.005x0.5 


L z 


p 

"to 


k D S D 


60x0,005x0.5 


l f 


0,1 


k- 5 r 


165x0.01x0.5 



= \ jn k/w 

= 1.333 K/W 
= 0,121 A/fK 



Finalement: 

o= 



1850-25 



0,111 + 1,036+0,121 



= 143' 



2.7.3 Cylindre creux a surfaces laterales isotne-r-es 




Ti 



T z 



j t i-h t « 



Figure 2.10: Cy "i3 r e e'e-.x 5 s 
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L'equation generate de ia conduction en coordonnees cylindriques s'ecrit (c.f. 
equation 2.15): 



dT ( S 2 T 1 dT \ d 2 T d 2 T) 

dl \6x r dr r SB' 8z' ) 




Dans le cas ou le transfert se fait dans une seule direction (r par exemple), ne 
depend pas du temps (didt = 0) et qu'il n'existe pas de sources de chaleur (£)' = 0), 
l’equation devient: 



d 2 T \ 8T 
d r 2 + r dr 



(2.30) 



avec !es conditions aux limites: 

Mr=ii)=7; 

[. T[r = cj = T. 



L'equation (2.30) peut aussi etre exprimee com me: 



d_ 

dr 




= 0 



et sa solution est: 



r 



dT 

dr 



= C, 



T=C } Ln{r) + C 2 



Les conditions aux limites: 



Jr = r, : T x = C 1 Zff(^) + C 
Jr = r 2 : T-C x Ln{r^) + C, 

permettent d’exprimer les constantes d’integ ration: 



C, 



T-T, 






Ln 






T-T, 



Ln 



Vr,, 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



La temperature sera decrite par [’expression: 

T = C l Ln(r)+C. 2 
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vient: 



En remplagant dans I’expression de la temperature decrite par (’equation (2.32), il 



T[r) = ~ f T \ Ln O’) + 7J - Tl r\ Ln ( r i ) 



Ln 






Ln 



\r 2 J 



qui peut s'exprimer sous la forme: 



T= T t + ~ l A Ln 

Lnl^ 

VrJ 



A 



(2.35) 



On peut deduire I'expression de la densite du flux themnique: 



dT T x -L 1 
<P = -k-r = ~k~ * 



dr 



Ln 



\LJ 



(2.36) 



ainsi que celle du flux thermique: 



(£>= (pS = -k 



T, -T 2 1 



Ln 






V 2 J 



— IttLt = InkL 



T-Z 



Ln 



a 



(2.37) 



et aussi celle de la resistance themnique pour un cylindre creux: 

Ln\ 



A 



R = 



InkL 



(2.38) 



Exemole 2.4: 

Une cheminee en bdton anmri (k } = 1,1 W/m°c) possede un diamStre intSheur 
d l = 600mm et un diamStre extSrieur d e - 1000mm doit Stre revalue de tints rieur par un 
matSriau rSfractaire (k 2 - 0,5 W i m°c). DStermlnez: (a) I'Spaisseur du gamissage, (b) la 
temperature de (a surface extSrieure de la cheminSe pour que les pertes themniques ne 
dSpassent pas 2000 W / m et que la temperature de la surface intSrieune de la pami en 
bSton armS ne dSpasse pas 200° C. La temperature de la surface interne du gamissage 
est prise Sgate £ 425° C. 
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Solution: 



T -T 

a) Le flux de chaleur a travers le garnissage s'exprime: 0 = --- — 

R 



Ln 



avec: 



R = 



d) 



2 nk , L 



*2, 7i h T 

_e flux par unite de longueur s’exprime done: Olm ? — - 





(d\T 


Ln 


_L p 




\d) 



et entraine que: 



</ = 



600 



(425-2<»;2x*-0.5 



2i»0 



= 421,36 mm 



aboutissant a: 



600-421,36 on 
e = 89 mm 



(b) Le flux etant conserve, la quantity de chaleur traversant la paroi en beton arme 
s'exprimera: 



Q = 



2nL 



Ln 



■(t -T ) 

| \ F W 



ce qui entraine que: 
Ln 



T-T 

V 



( d . A 



V«i J 



2nk, 



= 200 -- 



( 1000^1 

Ln x2000 

V 600 ) 

2x7rxl,l 



= 52°C 
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2.7.4 Sphere creuse & surfaces isothermes 




Figure 2.11: Sphere creuse & surfaces isothermes 



En coordonnees spheriques, I'equation generale de la conduction s'ecrit (cf. 
equation 2.16): 



dT 


\d 2 {rT) 1 1 3 


1 

F-i 

'Tj 

i 


1 d 2 T 


Q 


~dt" a \ 


r dr 1 T r sin(^) dy/ 


L MUW ^j 


j ■ 

r 1 sin 2 (^) dtp" 


+ pc 



Dans le cas stationnaire {dl 3t = o) , unidimensionnel [d I dy/ ~ d I d<p = o) et sans 
sources (Q' = o) , I'equation se reduit a: 



1 d 2 {r T) 
r dr 2 



(2.39) 



qui peut s'ecrire: 



r dr 



dr 



(r T) 



\_d_ 
r dr 



dT 

T+ '* 



1 ( dT dT d 2 T 

- — + — + r—~r 
r\.dr dr dr\ 



On aboutitfinalement a: 

d 2 T 2 dT n 
dr 2 r dr 

avec les conditions aux limites: 

'7tr = r,)=j; 

‘.7lr = r ; ) = r ! 

La solution peut etre recherchee en posant: 



(2.40) 



(2.41) 



L'equation se transformer alors: 
rU' + 2U= 0 



(2.42) 





2 Transfert de chaleur par conduction 



21 



La solution recherchee aura la forme: 




d’ou: 



/•(/•) = -— +C 

r 



Les conditions aux limites: 

r 

^ = - — + C, 



et 



T =-^+C., 
‘ r. 



permettront de trouver les expressions des constantes: 

Cl ~ ~ j_ J_ 

C = £ - 



t x -t 2 


i r 




V] ^ ^ 


r \ 



Texpression finale de la temperature sera: 

I 1 

7W = 7; +(j; - r,) -j — ^ 7 

r. r, 



Trouvons I'expression du flux thermique. II s’exprime : 

<rr 

c P=-kS — 

Jr 



«vec 



JT 

dr 



T\-Z 



- A 



1 

r. 



1 



Fnalement: 



<D=-£S’ 



, 7;-r, 



i i 



z' i ’N 



v /w 



r, 



(2.43) 

(2.44) 

(2.45) 

(2.46) 

(2.47) 

(2.48) 

(2.49) 
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Exemple 2.5 : 

Un ballon utilise en laboratoire et contenant une huile moteur est sou mis S une 
source de chaleur (voir figure ci-dessous). Le ballon en verre (k = 0,80 W i m"C) est 
a ssimile a une sphere dont les diamStres intSrieur et exterieur sont respectivement Sgaux £ 
20cm et 21cm. L'huile et fair environnant possedent des coefficients de convection 
respectivement Sgaux a 1 5 W l m 2 °C et 10 W i nr'C . 



■;t - Ballon 






Source de chaleur 



II est demands de calculerle flux thermlque que doit foumir la source de chaleur afin 
de garder la temperature de l'huile constante a 80 : C. On prendra la temperature a mbiante 
Sgale S 20° C et on nSgligera les pertes parle haut du ballon. 



Solution: 

L’apport de chaleur necessaire pour maintenir la temperature de l’huile constante a 
80° C doit etre au moins egal aux pertes subies par le ballon. 




Le flux s’exprime (c.f. figure ci-dessus): 



o = 



Th-Z 






avec; • 



K,2 = 



Finalement : 



1 


1 


\A.n.r{ ' 


0 

" 15.4. ;r.(0,l) 2 


1 


1 


h a A.n.r{ 


" 10,4.*- (0,1 05) : ' 






Hr, r 2 J 


1 


kA.fr.r 2 


" 0,8.4, *-.(0,l) : " 



= 0,530 °CfW 
T = 0,722 °C./W 



T = 0,520 °CIW 



tf> = 



80-20 



0,530 + 0,052 + 0,722 



= 46 Watts 
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2 8 Probleme de conduction unidimensionnelle avec d6perdition a 
travers les surfaces laterales 

2 8.1 Equation generate 

Ce cas represente le probleme type de conduction dans un solide (refroidissement 
ca' ailettes, fils chauds d’anemometrie etc.). II existe dans les corps longs ( x yy y etz) et 
e transfert se realise par la surface laterale. 



/ 

/ 

* r ^ 


* ^ 


. 1 . -ifJ. i x±dx 




/ 

/ 

/ 


1 1 


LJ 




1 





Figure 2.12: Bilan thermique a travers un element infinitesimal de tige mince 



Soit la tige representee en figure 2.12 dont la longueur est assumee tres grande 
oevant ses dimensions laterales. Elle possede une extremite chauffee en permanence. La 
cfialeur se propage dans la tige et se deverse dans le milieu ambiant. 

Sotent: 

• p : le perimetre de la section droite de la tige, 

• k \ sa conductivity thermique, 

• h: son coefficient d'echange convectif, 

• T a . la temperature ambiante, 

• S: la surface de la section droite de la tige. 

Dans les differents cas etudies precedemment, la formulation generale n'introduisait 
pas les conditions de surface et les effets des echanges peripheriques. Afin d’etre en 
mesure de les prendre en consideration, considerons un element de volume ( Sdx ) et 
‘aisons son bilan energetique. Le premier principe de la thermodynamique nous dicte; 

+Q'Sdx = U (2.50) 



Appliquons la relation de Fourier pour chaque flux: 

~ kS i%) +k ~ h pdx [ 7(x) - T °] + Q' Sdx = 0 



-kS\ 



df 

dx 



kS\ 



dt 



dx 



d (dt 



— +— — dx 



dx\dx 



- hpdx[T(x) -T a ] + Q'Sdx = 0 



En developpant, simplifiant et arrangeant I'expression precedente, on arrive a 1'equation 
d£crivant Ie transfert a travers la tige representee en figure 2.12: 

r.] = 0 



(2.51) 
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avec; 




(2.52) 



Remarques: 

1- C'est une equation importante, souvent appliquee pour les tubes et les fils possedant 
une section constante avec h et k constants. 

2- C'est une equation differentielle de deuxieme ordre a coefficients constants dont la 
solution generale est recherchee sous la forme 



T(x) -T } - A ch(mx) + H sh(mx) 


(2.53) 


% 

+ 

5 

r 

II 

1 


(2.54) 



3- Les conditions aux limites permettront de determiner les valeurs des constantes pour 
chaque cas particulier. 



2.8.2 Solution elementaire pour les ailettes unidimensionnelles 




L'ailette representee ci-dessus peut etre ass~~ ee a ~~e tge courte. Tout le flux 
n'est pas dissipe par sa surface lateral mas sa sector o r ote extreme est egalement 
traversee par un flux convectif. 

L'equation generale decrivant ce phenomene est ce e ~etee 2 51 

^--m 2 [T(x)-T]= 0 
ax' 

et sa solution generale est: 

T(x) -T, = A ch(mx) + B sh(mx) 

Les conditions aux limites imposent 

- A x^O : T(0) = T, => A = T - ^ 

- A x = L : le flux convectif est ega a- t'3“5~* s ca' cc n duction a travers la 

section droite 
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done: 



h\ T{L) - J!,] =h(T {l - 7) )ch(mL ) + hBsh(mL) 

En procedant a la difference entre les equations (2.57) et (2.58) et en simpiifiant, on 
aboutit a I'expression de B: 

\~hch(mL) - kmsh(mL)]{ L- t) 

B = ~ 1 1 — — — (2.59) 

hsh(mL) + bnch(mL) 



d'ou la solution: 



yw-r„=(r,-rj 



-- ch{mL) + sh{mL) 

ch(mx) - sh(mx) 

— sh(mL) + ch(mL) 
km 



(2.60) 



En posant G = — , en simpiifiant et en arrangeant, on aboutit a I’expression de 
km 

I'equation generate decrivant revolution de la temperature dans une ailette: 



T(x) - T a Gsh[m(L - x) + ch\m(L - x)] 
T 0 - T a ch{ml) + Gsh(mL) 



Le flux total dissipe est represente par la somme des flux dissipes par convection a 
travers !a section droite extreme et la surface laterale. I! s' exprime: 



<3> , = 



hS\T(,l)-T,\+\f>l{T{x)-T,\dx 



(2.62) 



Dans le present exemple ou le regime est permanent, le flux total dissipe est egal a 
celui transfere par conduction a travers la section de I'ailette: 



$ 



total 



. dt 
= — 
ax 



= AmS(r 0 -T a ) 



X=0 



sh(mL ) + Gch(mL ) 
ch{mL) + Gsh{mL) 



(2.63) 



Finalement apres arrangement, le flux aura pour expression: 



$ 

^ to Lai 



tmSfc-Z) 



th{mL ) + G 
1 + G th{mh ) 



(2.64) 



2.8.3 Efficacite d’une ailette 

La qualite d'une ailette est caracterisee par le rapport de sa performance effective a 
celle d’une ailette ideate de temperature uniforme egale a sa temperature de sa base. Une 
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telle ailette degagerait un flux de chaleur de: 
et son efficacite sera definie comme: 






s = 



ti el e cho^ 



rf) 

maximum e <hang cable 



(2.65) 



( 2 . 66 ) 



e = km 



s(t T \ th ^ mL > +G 1 
' ° a h + Gth{mL)hpl{r (l ~T a ) 



En multipliant numerateur et denominateur par ( mjhp ) et en simplifiant, on aboutit a: 



s - 



th(mL) + — 
km 

hL 

mL + — th(mL) 
k 



(2.67) 



Exemole 2 . 6 : 

Soit une ailette droite de section cirvulawe de diam&re 2cm et de longueur 10cm 
exposee £ un ecoulement convectif avec h = 25W ! m 2 K . II est demand 6 de comparer la 
distribution de la temperature dans cette ailette, ainsi que I'ecoulement de la chaleur et les 
efficacies pour trois materiaux diff&rents: (a) I'aluminium avec k = 200W/mK , (b) I'acier 
inoxydable avec k = 16 W ! mK ef (c) le verve avec k = \,2W ! mK. 

Solution: 

• La distribution de la temperature est calculee grace a t'equation (2.61): 

T{x) - T a Gsh[m(L - x)| + ch\m(L - x)] 

T 0 - T, ch (mL) + G sh (mL) 

avec: 

l§T_ f 25. ^ .0,02 1 5000 

m ~\kS~\k.7r.(0,0]y V k St km km 



d'oCi: 





m 


mL 


G 


Aluminium 


5,00 


0,500 


0,025 


Acier 


17,68 


1,768 


0,088 


Verre 


64,55 


6,455 


0,323 
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et la distribution de la temperature (a 77 A7J,) peut done etre calculee. Les resultats d'un tei 
calcul sont resumes dans le tableau ci-dessous. 





x(m) 


0,00 


0,02 


0,04 


0,06 


0,08 


0,10 


Aluminium 


1,00 


0,957 


0,924 


0,899 


0,88 


0,877 


Acier 


1,00 


0,719 


0,530 


0,407 


0,33 


0,306 


Verre 


1,00 


0,275 


0,076 


0,021 


0,00 


0,002 




• L'equation (2.64) permet de determiner le flux de chaleur dissipee par I'ailette (par souci 
de comparaison, on considerera que le flux dissipe par I'ailette en aluminium est egal a 
I'unite, 100%), tandis que l'equation (2,67) permet de calculer les efficacites. Les resultats 
sont resumes dans le tableau ci-dessous. 

Expression du flux dissipe: = kmS{T 0 - T J x + Gt} ^ mL) 



h 

ihimL) + 

A7?2 

Expression de I'efficacite: s = ^ — — 

mL +—th(mL) 

n 





Flux: du. (%) 


Efficacite: £ 


Aluminium 


100 


0,96 


Acier 


55,9 


0,54 


Verne 


16,1 


0,15 



On notera que I'ailette en verTe possede la variation la plus importante a la base 
(0 < r < 0,06) tandis que son flux est le moins important 














2 Transfert de chaleur par conduction 



28 



2.9 Conduction bi-dimensionnelle, stationnaire, sans sources 
Application de la m£thode de separation des variables 



Pour un corps isotrope et homogene, en regime permanent et sans sources de 
chaleur, ['equation generale de la conduction se reduit a: 



d^T_ d 2 T 

dx 1 + dy 1 




( 2 . 68 ) 



Figure 2.14: Domains illustrant la conduction thermique bi-dimensionnelle 



Essayons d'integrer cette equation pour un domaine rectangulaire en utilisant la 
methode de separation des variables. Les conditions aux limites sont: 



T(x,y) = \ 



0 

0 

0 



im 



si x = 0 et 0 < y < b 

si x = a et 0 < y < b 

si y = 0 et 0 < x < a 

si y = b et 0 < * < a 



Separons les variables: 

T(x,y) = X(x)J(y) 



c 2 T drX c'-T drY 

TT = Y ~Y~T « TT = x ~rT 
ox dx cy ay 



Remplafons dans I'equation (2.68): 



1 d 2 X 

JtlhF 



1 d 2 Y 

~ = Cste = -a 



Y dy* 



II est remarquable que les equations ci-dessus sont egales malgre le fait que 1'une 
depend uniquement de x et I’autre uniquement de v. Cect ne peut evidemment etre verifie 
que si les deux etaient egales a une constante qu'on a nomme a. Ces egalites nous 
menent done vers le systeme d'equations suivant: 



d 2 X 

+ ZX=0 a vec X(0) = X(a) = 0 
dx~ 

d 2 Y 



(2.71) 



dy 



y~XY=0 avec 7(0) = 0 et 7(A) = f(x) 
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L’equation en x est un systeme appele de Siurn-Liouville. Ses solutions non nulles 
sont recherchees sous la forme: 

A' - Ae"’ x => = mAe"' x et ~~ = m 1 Ae’ ,,s 

dx di r 

L'equation en x devient: 

m 2 Ae m +AAe mx = 0 

(m 2 +A)Ae mt = 0 (2,72) 

x appartenant au domaine [0,a] , e ne peut etre nulle et done: 

m 2 +k=0 (2.73) 

Trois cas se presentent: 

ler cas: k< 0 

On pose: k = -k 1 => m 2 -k 2 =0 => m=±k 

La solution aura la forme: 

X = A^ + A 1 e tx 

Recherchons les valeurs des constantes en appliquant les conditions aux limites: 

J X (0) = 0 = + A n => A 1 — — 

jx(a) = A/“ + A 2 er ka = A,{e ka - X ka ) = 0 

done: A l = 0 et A 2 = 0 => X = 0 

Cette solution ne peut etre acceptee. 

2eme cas: k- 0 

d 2 X 

— — = 0 => X = A,x + A - , 
cfx~ 

Recherchons les valeurs des constantes en appliquant les conditions aux limites: 

(X(0)=0=A 2 

[jf(far) - A^a - 0 => A, = 0 ^ 

Cette solution ne peut etre acceptee non plus. 
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36me cas: X >■ 0 

On pose: X- k 2 => m 2 + k 1 = 0 => m = ±ik 

La solution aura ia forme: 

X - A x cos(£r) + A 2 sin(-fcr) (2.74) 

Recherchons les valeurs des constantes en appliquant les conditions aux limites: 

f*< 0 )=°=4 (2 751 

[x(a) = A, sin(^a) 



Afin d'aboutir a une solution non triviale, A 2 doit imperativement £tre prise differente de 0. 
Done: 



sin(Aa)=0 => ka-nn et k - — , n=0,l,... 

a 

Ceci introduit une infinite de solutions non nulles de ia forme: 
A' (! (x) = A 2 x = A 2 sin^X ;t /r) 



(2.76) 



(2.77) 



Remarques: 

1- X R forme I’ensemble des valeurs propres du systeme et X„ ses vecteurs propres, 
2 - {X H } forme un ensemble orthogonal, e'est a dire que: 

] X m (x),X tI (x).dx = 0 si m*n 

0 



Pour les valeurs de ['equation possede comme solution: 



Uy) = A n e K ~' a -' 

(nit \ (nn 

r.W = C.ch{—y) + D.*{—y 



(2.78) 

(2.79) 



Recherchons les valeurs des constantes en appliquant les conditions aux limites: 



( 



nn 



r(o) = c„ = o => r n (y) = D*J*{—y) 



(2.80) 



II existe done une infinite de solutions telles que: 



( dTT ') , ( nn ^ 

T(x,y) = X{x)J{y) = E n sin[— . xj -^[-^y ),n = 0,1,... 



(2.81) 




2 Transfert de chaleur par conduction 



31 



II existe done suivant la valeur de n une infinite de solutions qui satisfont les trois 
conditions aux limites x = 0, x - a et y ~ 0. La solution generale est done une combinaison 
lineaire de ces solutions particulieres. Elle s'exprime: 

La constante E n est choisie de telle fagon que la derniere condition aux limites 
s'exprimant; 

^£ M sin^— xj sh\ ^-bj - f(x) soit verifiee. 



Exemole 2.7: 

Soit la plaque ractangulaire illustr£e ci-dessous. II est demand e de determiner la 
distribution de la temperature dans cette plaque en fonction de x et y. 




Solution: 

En proc6dant & la separation des variables et en rempla?ant dans I'equation 
generale de la conduction, on aboutit au systeme de Stum-Liouville: 



d-X . 

- + A 2 X = 0 



dx 

<fY 

dy 



-A 2 r = o 



dont la solution generale s'exprime: 
j X - A cos(/br) + B sin(Ax) 
| r = Ce*+De~* 



•ntrainant I'expression de la temperature: 

T(x,y) = X.V = [A cos(vbr) + 5sin(Ax)][Ce ,,:V + De ^'\ 
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Les conditions aux limites permettront d'evaluer les constantes: 

T y=b = 0 => [ A cos(a x) + B sin(x r)][C + D] = 0=> C = -D 
< 7^ 0 = 0 A{Ce* + De ' ) = AC{e Ay - e'*") = 0 => A = 0 
r v „, = 0 => [Bsia(AL)]c(e* - e~*') = 2BCsin(A L)sh(Ay) = 0 



La solution de cette derniere equation impose que sin(AL) soit nul c’est-a-dire que 

, nn 

A = — avec n- l,2,...d'ou I'existence d'une infinite de solutions dont la somme est ausst 
solution de cette equation, Cette derniere s'exprime: 



T(x>y) = y C„ sin 



nnx 



sh ■ 



nn y 



V L 



La derniere condition aux limites permet revaluation de la constante C„: 

tin 7 



^ nnx | 

X c „ sln \-r~ S/1 



rr - J 



V L 



T ■ f 

r - sm l.T) 



et impose que seule la constante C, est necessaire. Elle s'exprime: 

T 



C , = 



4? 



et la distribution de la temperature s'exprime finalement comme: 



n*,y) = T„ 



xy 

Sft — - 



nl 

L 



. ( xx 
sm — 

V L 



2.10 Conduction unidimensionnelle en regime variable sans source 
Application de la transformee de Laplace 

2.10.1 Generates sur la transform^ de Laplace 

La transformee de Laplace notee T(x,p) d’une fonction T(x,t) estdefinie par: 

n 

T(x,p) = \e pt T(x,l)dt 
0 

avec: • p: nombre reel ou complexe, 

• 7(x,p): Image ou transformee de la fonction T(xj). 
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2.10.2 Application de la transformee de Laplace a ['equation 
unidimensionnelle de la chaleur en regime variable 

L'equation generale de la conduction s’exprime (cf. equation 2.13): 

ST k (S 2 T S 2 T S 2 f] O' 

~T~ ~ I n -> + s 2 + -5 1 + 

cl pc\Sx~ dy t ) pc 

En unidimensionnel (3 1 Sy- S I Sz = 0) et sans sources (O'- 0), l'equation se 
simplifie pour s’ecrire: 

ST k S 2 T S 2 T 

o t pc Sx cx 



ou. 



S 2 T 1ST 
S x 2 a S t 



(2.85) 



Nous sommes done confrontes a une equation differentielle aux derivees partielles. 
Appliquons la transformee de Laplace a l'equation precedente (2.85): 



rfj 



i 



-pt 



o 




ST 

St 



St = 0 



Line integration par parties permet d’aboutir a: 

e" pt T(x,t)St - — 
a 

En remarquant que les deux integrales representent en fait la transformee de 
Laplace et en simplifiant, on obtient: 




[e-' > nx,t)l+p\e-''T(x i t)dt =0 



ckc a a 



( 2 . 86 ) 



avec: T n - T(xfi) 



La transformee de Laplace a done permis la simplification de l'equation de la chaleur 
gui est une equation differentielle aux derivees partielles en une equation differentielle aux 
derivees totales. Un changement de variables adequat permettra d'eliminer le second 
membre de l'equation: 



r = r-T t \ 



=> r 





dC-T 

Sx 2 



fequation (2.86) deviendra: 
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dont la solution est recherchee sous la forme: 

T\x,p) = Ae-* + Be* 

„ p 

avec: k= — 
a 



( 2 . 88 ) 



A partirde la solution transformee T*(x t p), la solution originate 7*(xj) = T(x,i)-T lt 
peut etre obtenue grace a la transformee inverse (des tables donnant cette transformee 
sont disponibles). 



Exemole 2.8: 



Une sphere en acier de 10cm de diamdtne et initialement £ la temperature uniforme 
de 300 a C est soudain plong4e dans un milieu o u rdgnent une temperature de 100 e C avec 
un coefficient de convection de \QW fm 2 °C. fl est demand 6 de determiner le temps 
nGcessaire a la sphere d’atteindre une temperature de 150 a C. 

CaractSristiaues de la sphere; 

p - 7850£g / m ' , c= 0,46 kJ l kg°C, k = 46 W i mK 



Solution: 

En assumant que la temperature de la sphere reste unifomne durant le processus de 
refroidissement et en negligeant les pertes de chaleur par conduction, celles echangees par 
convection seront manifestement representees par une diminution de I’energie interne de la 
sphere consideree. On peut done ecrire que: 

o = /.i(r-r„)=-pcV^ 



En assumant qu'au temps t=0, la temperature T = T 0 , il vient: 

T~t --^7. 

— = p > vV 

T - T 

*0 • £ tp 

avec: 7 = 150°C, £ = 100°C, T 0 = 30 CPC 



Done: 



hS 

p ct> 



10.4. ff. (0,005)- 
7850.460. (4/ 3). tt. (0,005) 3 



= 1 ,661.1 0^s~' 



d'ou: 



150 100 __ ^.1.6(51.10^,1 

300-100 " 



Finalement: t = 8346s = 2,32 Heures 
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2.1 1 Mythodes numeriques de resolution de liquation de la chaleur 
2.11.1 G6n£ra Iit6s 

Dans beaucoup de cas pratiques, la geometrie du corps ainsi que les conditions aux 
limltes sont trap complexes pour permettre une solution analytique du probleme. On a alors 
recours aux methodes numeriques. Ces demieres sont des methodes approximatives 
basees sur la technique appelee es differences finies. 

Afin d’appliquer une methode numerique quelconque a un probleme de transfert de 
chaleur, certaines mesures doivent etre prises, le but etant d’approximer I’equation 
differentielle et les conditions aux limites par un groupe d’equations algebriques. Ceci est 
gen6ralement accompli en subdivisant le domaine continu etudie en une serie d’elements 
finis. 



Dans chaque element, on assume que la temperature est constante et correspond a 
celle de son centre. On remplace le systeme par un reseau de tiges joignant les centres 
des diff6rents noeuds. 

Si on assume que la conductivity thermique de chaque tige est la meme que celle 
du materiau, le transfert thermique a travers ie reseau de tiges joignant les differents 
noeuds sera une approximation acceptable du systeme considere. 

Si N noeuds sont choisis, on aura N equations. Le systeme forme par ces demieres 
peut Stre r£solu en inversant Ea matrice de transfert ou en utilisant des methodes 
numeriques qui consistent essentieltement k estimer les temperatures aux noeuds et a les 
corriger jusqu'a I'aboutissement a une solution acceptable representee par I’egalite (a une 
certaine precision) des flux entrant et sortant du systeme. 



2.1 1.2 - Etablissement des Equations aux differences finies 
en regime permanent 

Le plan (x,y) est decoupe en un reseau de mailles. Chaque element du reseau est 
sous (’influence des 6lements qui I'entourent (c.f. figures 2.15 et 2.16). 




i-l 



i+1 




i.y+i 



H ) i +1 



Figure 2.15: Discretisation du domaine 
en elements finis 



Figure 2.16: Illustration des differents nceuds 

ayant une influence sur le noeud central 
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Considerons I'element P fJ appartenant au domaine et possedant une epaisseur 

unite. 



Ar 

H ■■■»! 




Figure 2.17: Bilan thermique du noeud central 



Les equations aux differences finies pour les nceuds situes a I'interieur du domaine 
sont obtenues en realisant un bilan thermique de I'element de volume de section (ArA>’) et 
d’epaisseur unite qui entoure le noeud P UJ : 

A>‘ + {<P^_, + <p !+ ij) Ar + Q'tefy = 0 (2.89) 

La loi de Fourier permet d'exprimer les differentes densites de flux de chaleur: 



■’-I 

[ 

j 


T. 


- T 

Aj -i 


ArA y 




Ar 




X. 


- T 


* <P,] -' Ar Ay 


■ ~ - k — 


Ar 




T i- 

. _ l Ki 




" AxA y 


— 


Av 




T . 

- t-LL 


- T 
Vu 


* ArA> ; 


— ~K 


Av 



Dans le cas particutier d’une maille carree (Ar- A_y - A/), le remplacement des 
differents flux par leurs expressions respectives dans 1'equation du bilan thermique et une 
manipulation adequate permettent d’obtenir 1'equation generate aux nceuds internes qui 
s’exprime: 



% 






J-l 



+ + T t . u + - 42;;, + = 0 



k 



(2.90) 
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Remarque import ante: 



De la meme maniere, on peut obtenir des equations generates aux nceuds externes 
ou coins (i.e. situes aux coins du domaine considere). Ceiles-ci s’expriment: 



Pour un noaud externe: 




+ 217 1 - 



h.AL 



T -—{iT. , + T 




( 2 . 91 ) 



Pour un Noeud coin: 



h.AL 



+ l\Z-2 









( 2 . 92 ) 



2.1 1.3 Conditions aux limites 

Trois conditions aux limites peuvent etre envisagees. 



2.11.3.1 Temperature impos^e 

Cette condition n’introduit pas d’equations supplemental. II suffit dans ce cas de 
reporter les valeurs des temperatures imposees aux noeuds correspondants appartenant a 
la frontiere dans !e systeme d 1 equations relatives aux nceuds internes. 



2.11.3.2 Density de flux imposge 

Cette conditions introduit une equation supplemental re. 



9 




Figure 2.18: Illustration de la condition limits de ta density de flux impos£e 
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Pour ('element de surface frontiere considere et si q> represente la densite de flux 
imposee a la frontiere, le bilan thermique s'ecrira: 



+ [<p + + Q'^y~ = 0 

En appliquant la loi de Fourier et en arrangeant, on obtient: 

A/ (A /) 2 

^ ^ +27^ -4?: ; + 2 = o 



(2.93) 



(2.94) 



2.11.3.3 - Condition mixte 

Dans ce cas, le flux impose et note <p precedemment sera lie a la temperature 
imposee qu'on notera par la loi de Newton qui s'exprime <p = - 7J ; ) . En remplafant 

q> par son expression de liquation precedente et en arrangeant, on aura: 

( AA A/ (A /) 3 

T 1J+l + T u _, + 2 T 2j - [4 + 2hj] T Ll +2hji;, + Q'^=0 ( 2 . 95 ) 



En 1'absence de sources de chaleur, cette relation peut etre simplifiee de la maniere 
suivante : 




+ T 



i.j-i 




= 0 



d'ou : 



T 




Bi+2 



avec: 




est le nombre de Biot. 



(2.96) 



2.1 1.4 - Methode de resolution- Etude d’un cas pratique 
2.11.4.1 -Gen6raiit6s 

La discussion precedente nous a montre qu’une methode numerique de resolution 
n’est en fait qu'un precede permettant une approximation de la distribution de la 
temperature qui est en realite continue par sa discretisation en un certain nombre de points 
appeles nceuds. II va de soi que plus le nombre de nceuds est important, plus on approche 
la solution exacte. Neanmoins, ceci implique aussi plus d’equations a resoudre et par 
consequent plus de complexity dans la solution du systeme. 




2 Transfert de chaieur par conduction 



39 



Dans la pratique, le choix d’un grand nombre de nceuds peut ne pas Atre necessaire 
et I’utitisation d’un maillage peu dense pourrait aboutir a la mfeme precision dans les 
resultats que celle d’un maillage contenant un grand nombre de nceuds. Ainsi, il est 
vivement conseille de ne pas fortement discretiser des le debut de la recherche de la 
solution mais plutot de choisir progressivement des valeurs pour A / de plus en plus petites 
et d'observer le comportement de la solution. Si ie probleme est correctement formule, les 
temperatures aux nceuds devraient converger vers la valeur exacte avec la diminution de 
A/ . Notons que les erreurs quand le nombre de noeuds augmente sont cumulatives, 

Afin d’introduire les methodes de resolution les plus representatives, interessons- 
nous au cas d’un reseau dont les temperatures aux noeuds frontieres sont assumees 
connues (figure 2.19), II est demande de determiner les temperatures correspond antes aux 
nceuds internes. 



1.4 e 



(20) 1,3 f 



( 20 ) 1,2 f 



(10) (10) 

2.4 3.4 



e < 


1 < 

2,3 


> © 

3.3 


U 1 


) i 

2,2 


> 0 

3.2 


O " i 

0 « 


> 

1 


) ( 

* A 



A 4,3 (40) 



A 4,2 (40) 



A 4.1 



2,1 

(30) 



3.1 

(30) 



Figure 2.19: Illustration d'un domaine discr6tis6 dont les temperatures 
aux noeuds frontieres sont connues 



En appliquant I'equation aux noeuds internes demontree auparavant, on obtient: 
Noeud (2.2): T 2) + T 2i + T l2 + T i2 - 4 T 22 = 0 
Noeud (3,2): T u + T 2i + T 22 + T 42 - 4T i2 = 0 
Noeud (2.2V. T 22 + T 1A + T u + T n - 4T 13 = 0 
Noeud (2.2): T i2 + T 34 + T 2i + T A1 - 4 T 21 = 0 

On aura le systeme: 

-4 T 22 + T 23 + T r _ =-(T 2[ +T 12 ) 

T 22 ~4T 32 + T 3i = — (zj, + £n) 

T 22 -4T 22 + T 3i =-[T M + T l3 ) 

T 2i +T i2 - 4T 23 = ~(t 2A + 4 ) 
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2.11.4.2 - Methode matricielle de resolution 

Le systeme precedent exprime sous forme matricielle donne: 



(4 


1 


1 


O'' 




(t 22 ) 




l '-5Q > 


1 


0 


-4 


1 








-70 


1 


-4 


0 


1 








; -30 


lo 


1 


1 


- 4 ) 




l TJ 




V- 50 J 



qui peut s ecrire: 

U][T} = [B] 
dont la solution est: 

[ b ] 

La difficulty essentielle pour cette methode reside dans I’inversion de la matrice [4] . 
En effet, 1‘importance du nombre de naeuds rend complexe la solution qui demanderait 
alors des temps et des moyens de calcul importants. Dans ce cas, cette methode n'est 
plus pratique. 



2.11.4.3 - Methode des residus 

Une methode pemnettant un calcul a la main quand le nombre de noeuds est assez 
r6duit est celle des residus. c’est une methode iterative qui consiste a approcher la solution 
(a la precision demandee) en procedant de la maniere suivante: 



lere 6tape . Se fixer des valeurs pour les temperatures inconnues T* ]) , 



2eme etape : Introduire ces valeurs dans le systeme d’equations et calculer le residu de 
chacune d’elles ^ a) , 

3eme etape ; Chercher £ eliminer les differents residus en commengant par I’equation qui 
possede la plus grande valeur de ce dernier en valeur absolue et en 
remplafant la valeur de la temperature inconnue qui lui correspond par: 



Tg 2) = + AT avec: AT = - 



pm 



coefficient de ] 



Conclusion : Poursuivre le calcul jusqu’d ce que tous les residus 

deviennent suffisamment petits pour n’entrainer que des 
modifications negiigeables des temperatures. 



Example 2.9: 

Re prendre le rdseau de nceuds repr£sent£ en figure 2.19 et procdder £ sa resolution 
par application de la methode des n&sidus. Une precision de 0,2“ est requise. 
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Solution: 

En appliquant I’equation g6n6ra!e aux differences finies aux ncBuds internes, on 
obtiendra le systeme suivant: 

4T 22 - T 2 ~ - T n - T 2] - T u = 0 
44-T 22 -r S3 -r 31 -4 = o 
'4T 23 -T 12 -r i ,.T 2A -T u = 0 
^-T 2i -T yi -r M -r« = o 

avec les conditions aux limites : 

T l2 = T u = 20 (T u = 25 

^4 = T » = 1 ° 1^4 = 15 

'T 42 = T a2 = 40 8 7^ = 25 

7 2l = ^ = 3o It;, = 35 

La resolution de ce systeme commence par le choix de valeurs pour les 
temperatures inconnues. Soit: 

7™ = 25 ; 2™ = 35 ; 7™ = 15 ; 7™ = 25 

Remplaqons les temperatures dans le systeme d'equations precedent par leurs 

valeurs respectives. On obtient: 

7?, (1) =0 , R™=0 ; R? ■ = -20 ; R™=0 

Corrigeons la temperature qui correspond a ('equation qui a eu le residu le plus 
important: 

of 1 ) on 

AT™ = --2- = — = -5 => £?> = 7l n + A 7™ = 35- 5 = 30 
4 4 

Les nouvelies valeurs des temperatures seront: 

7^ = 25 ; 7^ = 30 ; 7™ = 15 ; 7™ = 25 

Calculons les residus: R^ v = 5 ; R 3 2) = 0 ; R 3 2) = -20 ; Rj 2) = 5 

Corrigeons T 2% \ 

Rl 2) 20 

A 7™ = _i^L— = _ = 5 jo) = + AT™ = 15 + 5 = 20 

4 4 

Done: 7™ = 25 ; 7™ = 30 ; 7™ = 20 ; 7™ = 25 

Calculons encore les residus : 

7?™ = 0 ; R‘ 3) = 0 ; R™ = 0 ; R™ = 0 
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Les Ry calcules a partir de cette iteration seront nuts. La precision de 0,2° 

demandee est atteinte (dans ce cas particular, les valeurs des temperatures atteintes sont 
des valeurs exactes) et le resultat final est: 

T 22 = 25 ; T, 2 = 30 ; T 2i = 20 ; T i3 = 25 



2.11.4.4- Ntethode iterative de Gauss-Seidel- Utilisation du principe 
de resistance thermique 



Lorsque le nombre de noeuds est important, I'utilisation d'une methode numerique 
iterative devient necessaire. L'une de ces methodes est cede de Gauss-Seidel qui peut 
etre appliquee par I'utilisation du principe de resistances therm iques etudie a la section 2.7. 
En notant notre noeud 'i ’ et ceux I'entourant ’ / ' (c.f. figure 2.20) et en effectuant un bilan 
thermique, il vient: 






Tj-T, 

% 



= 0 



• q t : Chaleur delivree au noeud V ’ par une source, unrayonnement etc., 

• Ry : Resistance thermique dependant du mode de transfert entre les nceuds 'i ' et 

' 7 ' {conduction, convection aux frontieres etc.). 




Figure 2.20: Roseau de noeuds 'j' entourant un noeud interne Y 



Remarque: 

En remplaqant le second membre de cette equation par un residu, celle-ci peut etre 
resolue par (’utilisation de la methode de relaxation. 

La temperature au noeud interne ' i ’ s’exprime: 
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Cette equation est utilisee dans I’algorithme de Gauss-Seidel qui est applique selon 
le principe suivant: 

lere Etape: Se fixer des valeurs initiates pour fes temperatures T., 

2eme etape: Les nouvelles valeurs des temperatures sont calculees par application de 
Itequation (2,95) en utilisant toujours les dernieres valeurs des T., 

3eme etape: Le calcul se poursuit jusqu'a ce que la difference entre les valeurs des 
temperatures aux etapes ( n ) et (w-1) soit negligeable c'est-a-dire que 
\jin) _ j’fn+uj < s s etant un nombre assez petit. 



Exempt e 2.10: 



Reprenons le mime exemple traite pricidemment gr&ce aux mithodes matricielle 
et de relaxation . En nofanf les nceuds internes 1,2,3 et 4, celui-ci se presenters selon la 
figure 2.21. 




Figure 2.21: Illustration du r6seau de nceuds pour application 
de Gauss-Seidel 



Solution: 

Dans ce cas, le transfert thermique entre les noeuds se fait par conduction et aucune 
source de chaleur n'existe. En plus, le materiau est assunrte isotrope {k- Cste). Ceci 
entraine: 

Ay Ax AI 1 
^ k Ay kAx k A! k 
c h = 0 

doii: 




2 Transfer! de chaleur par conduction 



44 



Appliquons I'algorithme de Gauss-Seidel en choisissant: 

* « = 0 T x = 15 , t = 25 , T. = 25 , T> = 35 



7; = ^(20 + 10 + r, + r 3 ) = ^-(20 + 10 + 25 + 25) = 20 
% =^(lO + 40 + ?; + 7 4 )=^(lO + 40 + 20+35) = 26,25 

= ~(2o + 30 + 7; + t 4 ) = ^(20 + 30 + 20 + 35) = 26,25 
T a = ^(40 + 30 + r, + 7;) = ^-(40 + 30 + 26,25 + 26,25) = 30,625 

?; =-^-(20 + 10 + 26,25 + 26 , 25 ) = 20,625 
Z = “( 10 + 40 + 20,625 + 30 , 625 ) = 25,312 
T 3 = ^-(20 + 30 + 20,625 + 30 , 625 ) = 25 , 3 12 
T 4 =^(40 + 30 + 25,312 + 25 , 3 12 ) = 30,156 



Les calculs se poursuivent jusqu'a I'obtention du resultat a la precision souhaitee. 
Le tableau ci-dessous donne les resultats des quatre premieres iterations avec une 
precision inferieure a 0,12°, 



Iteration 




r. 


t 3 


T, 


6 mm. 


0 


15 


25 


25 


35 


- 


1 


20 


26,25 


26,25 


30,625 


5 


2 


20,625 


25,312 


25,312 


30,156 


0,938 


3 


20,156 


25,078 


25,078 


30,039 


0,469 


4 


20,039 


25,019 


25,019 


30,009 


0,117 



2.1 2- Regime variable 
2.1 2.1 - G6n6ralit6s 

En regime variable {d!dt~ 0), unidimensionnel {didy = dldz- 0) et sans 
sources internes, ['equation de la chaleur s'ecrit: 

dT k d 2 T 
dt pc dx 2 

ou: 

d l T 1 dZ 
dx 2 a dt 
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avec: 



a = — : Coefficient de diffusivite thermique. 

pc 

Le calcul de la temperature peut etre conduit de deux fagons; 

1- En exprimant la temperature T au temps ( t + At) en fonction des temperatures connues 
c'est-a-dire celles au temps t : Cette methode est dite explicite. 

2- En exprimant la temperature T au temps (/ + At) en fonction des temperatures 
inconnues au meme temps ( t + Ar) : c'est la methode implicite. 



2.12.2- Methods explicite 



Exprimons I'equation de la chaleur en regime variable en termes de differences 
finies. Pour cela, utilisons le developpement en serie de Taylor : 



T - T + 



( d f\ 



\dx) 



t , IdT 
Ti ~'~ Ti \dx) 



Ax + 



Ax + 



<? 2 rl(Ax ) 2 



dx* J 2 
d 2 T\ (Axf 



.dx-) 2 



£,+£i = 27 ? + 



f d 2 f 



(Ax) 2 



d'ou : 



d 2 T 
dx 1 ~ 



(At)’ 



(2.98) 



En remplagant la differentielle de T par son expression dans [’equation de la 
chaleur, celle-ci devient: 



(Ax) 2 "a At 



(2.99) 



En posant: M - 



(Axf 
a At 



et en arrangeant, on obtient: 



( 2 . 100 ) 
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Remarque: 



L’avantage de cette technique est que chaque equation consideree ne peut contenir 
qu‘une seuie inconnue. Par contre, n'importe quelle valeur de M ne saurait convenir. On 
montre que la condition de stability pour un noeud interne impose: M> 2 fois le nombre de 
variables geometriques. 



Le meme probleme de stabilite se pose pour les noeuds situes sur la frontiere. On 



montre aussi que la condition de stabilite dans ces cas impose: A/S2 




2.12.3- M4thode implicite 

Pour cette methode, I'equation est developpee en tenmes en termes d'inconnues de 
T au temps ( t + Af) . Reprenant ['expression prec6dente, I'equation de la chaleur s'ecrit: 



y yf + Ar . - Af 

A'4-l *?-l 




(A Jr) 2 



t t ; +s ’ + r 

a A t 



( 2 . 101 ) 



En developpant I'equation et en arrangeant, on obtient: 

(m +2)t; + a> - t‘;*’ - 1?_\* = mt; (2.102) 



2.12.4 - Methode g6n£rale 

D'une fagon generate, I'equation de la chaleur en regime variable peut etre exprimee 
comme une combinaison de deux termes: Tun calcule a partir des expressions au temps t , 
et I'autre a partir des expressions au temps (t + Ar) , Pour un noeud interne, elle s'ecrit: 




1 

a At 





(2.103) 



• a = 0 On aura la forme de I'approche implicite, 

• a = 1 On aura la forme de I’approche explicite, 

• a - 0,5 On aura la forme dite de 'Crank-Nicholson'. 



Exercices 



2 . 1 : 

D4crine brieve meat le mScanisme responsable de la conduction de la chaleur dans 
un milieu solide. 
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2 . 2 : 

Comment d£finissez-vous ta conductibilite thermique (ou coefficient de conduction) 
d’une substance ? 



2.3: 

Comment d£finissez-vous la resistance thermique ? Quelle est la relation qui la lie £ 
la resistance elect rique ? 



2.4: 

Quelle est I'importance de I’utilisation de I’analogie existant entre les grandeurs 
thermiques et eiectriques ? 



2.5: 

Expliquez brievement les hypotheses entrephses a fin d'arriver £ la loi de Newton 
donnant I' expression du flux convectifdes surfaces. 



2 . 6 : 

Decline bri£vement le mecanisme responsable du transfert de chaleur par 
convection et d£duire ses principales differences avec celui realise par conduction ? 



2.7: 

Enum£rez les etapes essentielles de la mise en place d'une methode numerique de 
resolution de T equation de la chaleur. 



2 . 8 : 

En quoi rayonnement thermique differe-t-il des autres types de rayonnement faisant 
partie du spectre du rayonnement 6lectromagnetique ? 



2.9: 

Qu’est ce que I’emissivite d’une substance ? 



2 . 10 : 

Quelle est la difference entre un corps noir, un corps gris et un corps reel ? 



2.11 : 

Les deux faces d'une plaque de cuivre de 3 cm d'epaisseur sont maintenues a ux 
temperatures de 400°C et 100°C. Calculez la quantite de chaleur transmise par unite de 
surface a travers la plaque (k cuWn = 310W t m°C). 

R6p.: Q/ S = 3JMW / m 2 
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2 . 12 : 







1 Sou 


rce 1 



Wattmetre 




Isolation 

Echantillons 

Cliauffage 

Isolation 



Aftn de mesurer la conductibilite thermique, deux Echantillons de matEriau identiques 
de dimensions (15x15 cm 2 ) et d'Epaisseur 2,5 cm sont placEs dans le banc 
d'expErimentation neprEsentE ci-dessus. Du courant Electrique est foumi au systEme de 
chauffage. Un wattmetre, placE dans le circuit, indique que la puissance dissipEe est de 10 
Watts. Les thermocouples attachEs aux surfaces chaudes et froides des Echantillons 
indiquent des valeurs pour les tempEratures de 300K et 260K raspectivement. Calculez la 
conductibilitE thermique k du matEriau. 

REp.: k = 0,\4IV / mK 

2 . 13 : 

De I'air E 20°C souffle sur une plaque de dimensions (50x75 cm 2 ) maintenue E une 

0 

tempErature de 25Q°C. Le coefficient d'Echange par convection h est Egal E 25W im 2 C. 
Calculez la chaleur transmise par la plaque. 

R6p.: Q = 2,\56 kW 



2 . 14 : 

Un courant Electrique passe dans un fil de 1mm de diamEtre et de 10cm de 
longueur. Le fil est immergE dans de I’eau E la pression atmospherique et I’intensitE du 
courant est augmentEe jusqu’E ce que I’eau arrive E Ebullition. Dans ces conditions, le 
coefficient de transfert convectif est trouvE Egal E 5000W/m 2 °C. Quelle est la puissance 
Electrique qui doit Etre foumie au fil afin de maintenirsa surface E la tempErature de 1 14°C? 

R6p,: Q = Puissance = 21,98 Watts 



2 . 15 : 

Deux plaques infiniment longues et parallEles dont les conditions de surface 
approchent celles du corps noir sont maintenues aux tempEratures de 1070K et 520K. 
DEterminez le transfert de chaleur par rayonnement par unitE de surface entre les deux 
plaques. 
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2 . 16 : 



Flux 



T- 

i 



Briques refractaires 



it 



+- 



Briques isolantes 



Le mur d'un four est compose de deux couches. La premiere, en briques 
refractaires, possdde les caracteristiques L x = 0,20m et k x = \,3%W / m°C . La seconde 
couche, en briques isolantes , est caracterisbe par L, = 0,10m et k 2 - 0,\1W ! mC . A 
I'interieur du four, la temperature T t estde 1650°C. Le coefficient d'echange convectif sur la 
paroi interne est h t = 10W/m 2 "C. L'exterieur est constitue d'air dont la temperature et le 
coefficient d'echange convectif sont respectivement T c = 25° C et h c = 10 W I m 2 "C . Calculez: 

1- les pertes de chaleur par m 2 de surface, 

2- les temperatures de la face interieure, de I' interface et de la face exterieure du mur, 

3- les pentes des droites T(x) pour chaque partie du mur. 



R6p.: 1- $= 1917^ / in* 



2- int _ = 1622°C ; T mxeiftlct = 1344°C ; T faccm = 216 T 



3- 



df) 



Re fraclaics 



= -1390°C/m ; (f) 



= -11280°C/ m 



}xol<yii/cs 



2 . 17 : 

Utilisant I'analogie existant entre les grandeurs thermiques et eiectriques et 
assumant une conduction unidimensionnelle, calculez le flux de chaleur par unite de 
surface & travers le mur compose illustre a la page suivante. 

A.N.: k A = MSW ImK , k a = 35W!mK , k c = S0W f mK k D = 55W ImK 




' unite tie surfoce 



= 44 kW 
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2 . 18 : 

Soit un fil de 1mm de diam&ne couvert d'une couche isolante de 2mm d’Spaisseur et 
de coefficient de conduction k - 0,5W / mK. La temperature ambiante de fair entourant le 
fil est de 25°C ef son coefficient de convection h = WirrrK. La temperature de ta surface 
du fil est de 100 e C. On demande de calculer le flux de chaleur degag£ par le fil par unite 
de longueur: 

1 - sans couche isolante, 

2- en presence de la couche isolante. 

R6p. : 1 - Transfert par convection uniquement: <X> = 2,36 W / m 
2- Transfert par conduction et convection: = 10,9 W Im 



2 . 19 : 

Une conduite cytindrique isotee est reptesentee ci-dessus. A I'interieur de la 
conduite s'Gcoule un fluide chaud (T t ,h r ). L'isolation est reptesentee par un cylindre creux 
de rayons R. et It. On assumera un regime permanent et une conduite longue 
(L >^R x ,ILet R,). 




On demande de dSterminer: 

1- Le coefficient d'Gchange global par unite de longueur de la conduite , 

2- le flux de chaleur transmis de I’interieur vers I'extGrieur par unite de longueur de la 
conduite. 

Application num£riaue: 

i?, = 60mm, R, = 65mm, R. = 80mm, = 46 W l mK, k 2 - 0,\6W i mK 
ht = mW/rn 2 K,h e = 20W/m 2 K,Ti= 333 K,T t = 293£ 

R6p.: K tUM =5,9SW/m 2 K ; 2-®=120P*7m 

2 . 20 : 

Determine! les dGperditions thermiques O se tenant £ fra vers une surface vittee de 
surface lm J ainsi que la repartition des temperatures dans les deux cas suivants : 
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1- Vitrage simple d'4paisseur e = 4 mm, 

2- Double vitrage compost de deux lames de verve d'4paisseur 4 mm chacune et d'une 
lame d'air IntermSdiaire stationnaire d'4paisseur 6 mm. 

Application Numerigue: k verre = \,2fV i mK , k ail ~ 0,024W I mK 

h = \2W 1 rtrK ; 7? = 20°C ; T t = CPC 

R6p.: 1- 3> = 117,6^ ; Temps: 10,2/9,8° C 

2- 0> = 47,2 FT ; Temps: ]6,1/15,9/4,1/3,9°C 



2 . 21 : 

Un distillateur solaire (voir figure ci-dessus) dans tequel s'ecoule de I'eau 4 la 
temperature de 360K est expos4 au rayonnement solaire dont le flux thermique 4 t ravers la 
plaque de verve est de 490 W/m 2 . La temperature de I'airambiant est de 300K. 

En n4giigeant I'effet de ia plaque de verve sur ie transfert thermique et en 
considdrant que ie bas et les cdt4s du distillateur sont adiabatiques, dSterminez le temps 
necessaire pour transferer 10 h Joules 4 1m 2 de surface de I'eau. 



AIR 



X ^XX 7 XX a 



Rayonnement incident 

Plaque de verre 

Eau en contact avec 
la plaque de verre 



✓ 



a. 



Onprendra: />„ re = 28,3 Wlm 2 K ; h m ^ atr = 6,8 WfnrK 
Rfep. : t - \h42mn 



2.22 : 




Consid4rons une bane cylindhque de longueur L dont une extr4mit4 est chauffSe 4 
une temperature uniforme T a . L'autre extr4mit4 est consid4r4e adiabatique. Soierrt p le 
p4rim4tre de la section droite, S I’aire de cette section, k la conductibilit4 thermique du 
mat4riau, h le coefficient ailette-air et T a la temperature ambiante. II est demand4 de 
trouver la distribution de temp4rature ainsi que le flux total de dissipation. 

= JThpS - T^th(mL) 



Bit; nx)-T, = (T e -z) 



ch(mx) - 



sh (mL) sh(mx) 
ch(mL) 




2 Transfert de chaieur par conduction 



52 



2.23 : 

Trouvez la distribution de temperature T(x,y ) dans la longue bane rectangulaire 
dont la section est representee ci-dessous. 




Une tdle metallique d'epaisseur 2L, de longueur infinie et de temperature initiale T, 
uniforme se trouve piongde brusquement dans un bain isotherme e temperature constante 
T a . En admettant que les deux faces (x = -L et x = +L) prennent instantandment la 
temperature du bain (i.e. au temps t: T(Lj) = T(-L,t)= T a ), trouvez la distribution de 
temperature T{xJ). 

{In + 1) nx 

2 T~ 



R6P. : T(x,t) = 



41 



I 



IZIM-I) fr 

, 4 /- 



(2n + l)x„ = 



cos 



2 . 25 : 

Un profile en 1 de hauteur 300mm dont le dessus est maintenu e une temperature 
de 250°C et le dessous & 90°C possede une epaisseur e-12mm. . De I'air e 250 e C souffle 
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le long de ce profile entrainant une valeur pourle coefficient de transfer! h~40W/m 2 K. . La 
conductivity thermique de I'acier est k=45W/mK. 

1- Trouvez la distribution de la temperature le long de la tranche centrale du profile, 

2- Dessinez la courbe qui re prd sente cette distribution. 



► f«-]2min 



300min 



90X 



*'\25()“C . - 
12tmn 



300inm 

T 



X t 



Rep. : 1- T(x) = 250 -lj.sh(\2,4x) 
2- 



X 


0 


0,01 


0,04 


0,1 


0,3 


0,336 


0,4 


0,5 


0,6 


T ( x ) 


250 


249 


246 


238 


91 


0 


-299 


-1647 


-6305 




2.26 : 



<P, 



<P, 



(1) (2) 



0.5 



►K- 



l 



(3) 



Figure 1 



(4) (5) 



1 



. 0.5 , 



t 




(i-1) 


t 




n 


Jl 



<Pi 

Figure 2 



Une aube de turbine de longueur 4cm et de diametre 1cm est constitute d'acier 
(k=25W/mK). la temperature de sa base Tj = 750K. L'aube est exposte £ de i'air chaud 

(T, - i ^ = W4QK, h = 450 W / m 2 K). II est demande : 

1- d'etablir I'equation generate aux noeuds internes (notts 2, 3 et 4), 

2- d’dcrire les trois equations aux nceuds internes 2, 3 et 4, 

3- de d terminer le flux de chaleur (dissipd ou gagne) par l'aube. 
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Rep.: 1 - AT + BT , - BT + C = 0 avec: A = -hpDAx , B = — — — , C. = hpD tsx 
' 1 4Dx 

2- Noeudl. AT 2 + BT x - B T h + C = 0 , Noeudy AT, +BT 2 - BT 4 + C = 0 
Noevd.4 ; AT a + BT 2 ~ BT^ + C = 0 

3- O = 65,45^ 

2 . 27 : 



300K 




12 3 4 



Une cheminee (k=0, 7W/mK) possede une section droite telle que definie par la 
figure. Si la temperature de la face interne est de 520K et celle de la face exfeme 300K, 
calcutez en utilisant une methode num&rique la distribution de temperature (Precision de 
IK). 

R6d.: T 21 = 362,5 * ; T 22 = 331,25* ; T n = 362,5 * 

2 . 28 : 



1 2 3 4 5 6 




Calculez grSce a une methode numerique: 

1 - les temperatures aux 12 points equidistants illustrGs sur la figure, 

2- le fiuxde chaleur par m£tre d'epaisseur. 

Donnies: k = 1,7 W ! mK ; Fa/re usage de la sym&rie. 

Rfep.: 1- T 25 = 266 * et T, s = = 277* ; 2-0 = 1,05 k-W / m 
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adiabatique 



Le soiide bi-dimensionnel illustr& ci-dessus (k = 20W / m c C) possede une source 

interne g6n£rant une quantite de chaleur de 9QMW ! m\ !l est demand £ d'exprimer les 
Equations aux nceuds internes et extemes ainsi que celui situ6 au coin (On assumera 
T 4 = 20 °C). 

R6p.: T, 0 + T, + T 2 ~4T & =-470 ; T 4 + T 6 -AT, = -570 

T ia + T 6 +T 9 + T u - 4T Ki =-450 ; ^ + T 7 + T h) -4T U = -550 
2% -T 6 - 0,5 T 9 = -440 ; 2 % - - 0,57; - 0,57^ = -450 

27; 4 -T l0 - 0,5T U -0,5T l} = -450 ; 27;, - T n - 0,5T lA = -400 
2T n -2T u -T 9 =-450 

2 . 30 : 



40° C 



40°C 



40°C 



Soit le reseau de nceuds illustre ci-dessus appartenant £ un soiide de conductibilife 
thermique k = 10W / m°C . Sa panoi de dnoite est soumise £ un flux de convection. 

li est demand 6 de calculer les temperatures aux nceuds (2,2), (3,2) et (3,3). 

FteP. : Noeud interne (2,2): T 2l + T i2 + T 2i + T u - 4T 22 = 0 

Noeud externe (3,2): 2,3 T i2 - 0,3 T w - T 12 - 0,51, - 0,5I ( = 0 
Noeud coin ((3,3): 2,67^0,67; - 2T 2i - 7J, = 0 
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d’ou: 

- 4 T 12 + T r _ = ~{T n + T 2i + T 23 ) = -90 
< 2 , 37 ^- 7 ^- 0,5733 = 5 
2,67;, - 7 J2 = 80 

La resolution de ces trois equations entraTne: 

7, = 23, ]°C , T n = 28,3° C et ?; 3 = 39,6° O’ 



Chapitre 3 : 

TRANSFERT DE CHALEUR PAR CONVECTION 



• Introduction 

• Equations generales de I’hydrodynamique 

• Concept de la couche limite, ses hypotheses et ses equations 

• Donnees empiriques de la convection 

• Methodes devaluation de h 

• Analyse dimensionnelle 

• Formufes empiriques couramment utilisees 



3.1 Introduction 

La convection est un mode de transfert de chaleur qui se produit le plus souvent 
entre un fluide en mouvemeni et une paroi solide. Alors que la conduction etudiee au 
chapitre precedent peut etre consideree comme un transfert d’energie du a des 
mouvements microscopiques, le phenomene de convection est un transfert du a des 
mouvements macroscopiques. On distingue deux types de convection: 

• La convection libre ou natureile ou le mouvement du fluide est cause par un 
champs de forces interieur (gravite, gradient de densite, gradient de temperature etc.), 

• La convection forcGe ou le fluide est mis en mouvement par Taction d'un champs 
de forces externe (pompe, ventilateur etc.). 

D’une fafon simplifiee, Tobjectif principal de I'etude du phenomene de la convection 
consiste essentiellement a: 

• developper des methodes permettant revaluation du parametre h , coefficient de 
convection; 

• etudier les differentes formules empiriques utilisees. 



3.2 Equations generales de I'hydrodynamique 

3.2.1 Generates - Definitions 

3. 2.1.1 Milieu continu 

Un milieu materiel est dit continu lorsque toutes ses proprietes varient continument 
dans Tespace et dans le temps. En d'autres termes, les distances considerees sont 
largement superieures aux distances intemnoleculaires. 
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3.2.1. 2 Milieu homog^ne 

Un milieu materiel est dit homogene iorsque toutes ses proprietes sont constantes 
dans tout le domaine considere (par exemple, la masse volumique dans les milieux 
incompressibles). 

3.2.1. 3 Milieu isotrope 

Un milieu materiel est dit isotrope Iorsque toutes ses proprietes sont identiques quel 
que soit I'orientation. 

3.2.1. 4 Fluide 

Un fluide se definit comme un corps presentant une vitesse de deformation non 
nulle si on lui applique des contraintes tangentielles aussi faibles soient elles. 

3.2.1. 5 Fluide newtonien 

Un fluide newtonien constitue un milieu materiel continu, homogene et isotrope. 



3.2.2 Equation de conservation *’ 

3.2.2. 1 Equation de conservation de la masse ou Equation 
de continuity 



Le principe de conservation de la masse d'un volume V nous dicte I'egalite entre la 
variation de masse dans ce volume et le 'flux de masse' a travers sa surface: 

= -JJ p V n dS = - jff div { p V\dV 

Y Ot s V 



~^div{p V) =0 



(3.1) 



Pour un ecoulement stationnaire (Bldt = 0) et incompressible (p-Cste), 
['equation de continuity se reduit a: 

div V ~ 0 



dn &v dw 
dx cy dz 



(3.2) 



Pour plus de details concemant les equations fondamentales ou equations de 
conservation, consulter un ouvrage specialise. L'auteur conseille celui classique de H.W. 
Liepmann et A. Roshko intitule 'Elements of Gasdynamics' et publie par John Wiley et sons. 
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3.2.2. 2 Equations de conservation de la quantity de mouvement 
ou Equations de Navier-Stokes 

En appliquant le principe de conservation des moments rapporte au systeme de 
coordonnees cartesiennes, les equations de Navier-Stokes s’ecrivent: 



Du 

P 17, =F '- & 



8P f d 2 u d 2 u d 2 u 



+ M 



Dv dP 

Dw 



+ 8y 2 + dz 1 . 
r d 2 v 8 l v d'"v 



\dx' 



= F - 

Dt * dz 



dy 

dP ( d 2 



— + ■ 



— + 



dy" dz ~ 



iv d~ w d ‘ h- 



+ Pl 



yd x- 



dy- 






8z " 



(3.3) 



avec: 



D e 8 8 8 

• — = - — + u — + v — + w — — : Derivee particuiaire, 

Dt dt dx dy dz 

• p \ Densite, 

• F x , F F z : Forces exterieures, 

• P\ Pression, 

• p: Viscosite dynamique, 

• u, v, w: Composantes du vecteur vitesse. 



3.2.2. 3 Equation de conservation de I’energie 

En appliquant le principe de conservation de I'energie thermique deduit du premier 
principe de la thermodynamique, I'equation de conservation de I’energie exprimee 
relativement a un systeme de coordonnees cartesiennes s’exprime; 



DT 

" c ’-5t = k 



d 2 T 8-T d 2 T ) 
8 x 2 + 6 y 2 + 6 z 1 ) 



En termes de la fonction enthalpie, elle s'exprime: 



V{ 



V 2 



H, +— + Q = H, +-r~ + W 
1 2 2 



ou: 



• H : Enthalpie, 

• Q: Chaleur regue, 

• W: Travail net produit. 
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De fair s'ecoule d’un reservoir a t ravers une conduite de section rectanguiaire tel 
que re pr^ sente ci-dessus. Si la hauteur du tube de sortie Y esf considerGe petite 
relativement a sa longueur L, le profit de la vitesse £ la section de sortie s'exprimera par: 

*'( y) = ^r{y Y ~y~) 

En assumant le debit assez faible afin que la density de fair dans le reservoir et la 
conduite de sortie soit consideree uniforme \p * /j(jr ,_y,z)], determine! (a) le ddbit 
instantane a la sortie et (b) le taux de variation de la densite dans le reservoir en fonction du 
temps. 

Donndes : 

V = 0,283 m 3 ; Y = 0,00254m ; L = 0,3048m ; P mT = 1,379.10 5 N I nr 
T tml = 277,SK- U tl „, . = 30,48m fs- ^ = 2S7,Q4J ttgK 

Solution 

a/ Le debit instantane- du fluide est fonction de la vitesse et de I’aire de la section de sortie. 
II s'exprime done: 

= II P udS = II pjr{y Y - y 2 H ,dz =pjt$S (y y - 

‘ } sz>rti# ^sortie ^ 

2 pULY 

Apres integration et arrangement: m vo , fic = — — 



b / En considerant la surface de controie comme etant celle enveloppant le volume V 
compose du reservoir et de la conduite et en utilisant fequation de continuity, il vient: 



fJJ^ rfv+ /W= 0= 




c ?r dp 2pULY 

La density ne dependant pas de jt, y et z, done: — V = -m ortie => — = -■ — ~ — 

(7 1 (/’I J V 

dp , 

L'application numerique donne: m wf!e = 0,027 kg ! s et — - = -0,096£g / m s 
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Exemofe 3.2 

De i'aira la temperature de 300 a C est detendu isentropiquement jusqu'a atteindre la 
vitesse de 250m/s. Calculez la temperature de I'air £ cette vitesse. 



Solution 



L'utilisation de I'equation de conservation de I'energie (3.5) permet d’ecrire: 

/ 1 U 2 
H,-H, =C t (T,-T 2 )=-^- 



d’ou: 



T 

1 2 




300- 



(250) ; 

2 x 1006 



269 °C 



3.2.3 Le concept de la couche limite et ses hypotheses 1 " 



3.2.3. 1 Concept de la couche limite 

Une solution complete des equations de Navier-Stokes representant un ecoulement 
visqueux pose des difficultes mathematiques considerables. Une simplification de ces 
equations generales est done necessaire, En 1905, Ludwig Prandtl remarqua que pour la 
plupart des ecoulements, I'influence de la viscosite n'etait significative que dans une region 
tres mince situee au voisinage de la paroi solide. Le reste du fluide peut done etre 
considere comme non visqueux et done ideal. Cette region est appelee: la couche limite. 

Considerons un fluide s'ecoulant aupres d'une paroi solide (cf. figure ci-dessous): 



u.. 



X, 


















' Paroi solide 

Figure 3.1 : Couche limite au voisinage d'une paroi solide 

Si le fluide est newtonien, la contrainte tangentielle s'exprimera grace a la loi derivee 
par Newton qui peut s'ecrire dans ce cas: 



*= fi- 



ll* 



(3.6) 



ou: represente la vitesse a la frontiere de la couche limite. 



Pour plus de details concemant I'etude de fa couche limite, consulter les ouvrages 
specialises. L’auteur conseille vivement celui ciassique de H. Schlichting intitule; 'Boundary 
layer theory' publie par Mc-Graw Hill. 
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A I'interieur de la couche limite, U„ est liee a la vitesse U par la relation: 

ch<_ = U* 

Gy 5 

La contrainte tangentielle aura done pour expression: 

du 

T = jj — 

dy 

On d6finit aussi la viscosite cinematique comme: 

_ /f. 

P 



(3.7) 



(3.8) 



Exemp/e 3.3 

Deux plaques horizontales sont placees £ 1,25 cm I'une de I'autre. L’espace compris 
entne elles est nempli d’une huile de viscosite dynamique p ~ 1,4 Ns/m 2 . II est demande de 
calculer la contrainte tangentielle exercee par I'huile si la plaque superieure se ddplace § la 
vitesse de 2,5 m/s. 



Solution 

L'application de I'equation (3.7) aboutit a: 



du 2,5 



280 N/m 1 



3.2.3.2 Hypotheses de la couche limite 

En tenant compte des particularites propres a I'ecoulement d’un fluide au voisinage 
d'une paroi solide, des hypotheses applicables uniquement a I'interieur de la couche limite 
et permettant la simplification des equations du mouvement sont enoncees. Elles sont 
connues sous le nom d'hypotheses de la couche limite. 

1- Le fluide immediatement en contact avec la paroi est immobile: e'est la condition de non 
glissement qui pour un ecoulement bi-dimensionnel se traduit par: 

j?/(r,0) = 0 
jv(x,0) = 0 



2- L'epaisseur de la couche limite notee 6 et definie comme la distance a laquelle 
u - 0,99(7,, est consideree petite relativement aux autres dimensions caracteristiques de 
I'ecoulement, Neanmoins, cette epaisseur augmente quand on se deplace dans le sens 
de I'ecoulement (cas de la plaque de I'exemple 3.4). 
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3- A i'interieur de !a couche limite, la vitesse axiale (dans le sens de 1’ecoulement) notee w 
est tres grande devant la vitesse radiale notee v c'est-a-dire que lorsque y^S aiors 
u >-> v. 



4- A I'interieur de la couche limite, les variations de la vitesse dans !e sens de recoupment 
relativement a y sont importantes. Ce qui peut se traduire par: 



Gu Gu Gv Gv 

=> — — 
ay Gx Gx Gy 



(3.10) 



Ces differentes hypotheses valables uniquement a I'interieur de la couche limite 
simplifient sensiblement les equations du mouvement. II en resulte en particulier que pour 
un couche limite bi-dimensionnelle: 



GP GP dP 

Gy Gx dx 



(3.11) 



Pour un ecoulement bi-dimensionnel stationnaire et laminaire (cf. Paragraphe 3.2.4), 
ces equations s'§crivent: 



Gv Gv 
Gx Gy 



Gv Gu 
w— + v— 
Gx ay 



1 dP G 2 u 

' ~ + v - i 

p dx Gy" 



(3.12) 



avec: 

j ~ v — 0 ay - 0 

\it = U B (x) = U„(x) a y-8 



L'epaisseur de la couche limite dans ce cas s'exprime en fonction du nombre 
adimensionnel de Reynolds (c.f. paragraphe 3.2.4): 

5x 

5, =~r= (3.13) 

lamwmrc ' • 



• 5, : Epaisseur de la couche limite laminaire, 

• x. Distance de I'origine (debut du corps) au point considere, 

• Re x .: Nombre de Reynolds calcule sur la base de la distance x. 



3.2.4 Regime laminaire et regime turbulent 

En 1883, Osborne Reynolds conduisit une serie d'experiences classiques impliquant 
des ecoulements a I'interieur des conduites. Ces experimentations d6montrerent 
I’existence de deux sortes de regime distincts: un regime dit laminaire et un autre dit 
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turbulent. Reynolds reussit a determiner le critere d'instabilite qui gouverne la transition 
entre ces deux regimes. 




IV 



Tube 1 



K 



Tube 2 









Figure 3.2: Representation de [’experience de Reynolds 



Ceci a ete realise grace aTinjection d'un mince filet de colorant a I'interieur d'un tube 
d'essai situe dans un reservoir rempli d'un liquide possedant la meme masse volumique 
que le colorant (cf. figure 3.2), Quand le debit est faible, le filet colore reste etroit et 
parallele aux lignes de courant dans le tube: c'est le regime laminaire (cf. Tube 1 de la 
figure 3.2). Lorsque le debit est augmente au-dessus d’une certaine valeur critique, le filet 
colore commence a onduler puis tres rapidement, on observe comme un eclatement de ce 
filet qui semble alors occuper tout le tube: c'est le regime turbulent (cf. Tube 2 de la figure 
3.2). Reynolds appliqua une analyse dimensionnelle aux ecoulements en conduite et 
conciut que la transition prend place pour une valeur fixee d'un certain parametre pouvant 
etre interprets comme le rapport des forces d'inertie a celles de viscosite. En son honneur, 
ce parametre est depuis appele: le nombre de Reynolds. La signification physique de celui- 
ci peut etre demontree comme suit: 




Figure 3.3: Representation d'un element de volume fluide 



Force typiqve d'inertie 



Masse x Acceleration 



U 2 , . 

T - plru 



= p If x 
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Force typique de viscosite = Contra int e de vis cosite x Surface 

U 

- /( - x L“ = p UL 



Done: 



Re = 



Forces d' inertie 
Forces de viscosite 



pUflf 
p UL 



pUL 

M 



(3.14) 



Utilisant la viscosite cinematique notee et definie par: v 
Reynolds s'exprime: 



— , le nombre de 

P 



Re 



UL 

v 



(3.15) 



II est generalement admis que la valeur du Reynolds transitoire definissant le 
passage de la nature iaminaire de recoupment a celle turbulente est egale a 1 *': 

- Re,. = 2300 pour les ecoulements en conduite, 

- Re t = 500000 pour les ecoulements sur plaques planes. 



Exemple 3.4 

Une plaque plane lisse de forme carree de 2cm de cote est tenue immergee dans 
de I’eau de viscosite cinematique v ~ 10 6 m 2 ! s s'ecoulant a la vitesse de U = 30 cm! s. II 
est demande de determiner la nature de I'ecoulement d une distance de 50 cm du bordde la 
plaque ainsi que I'Spaisseur de la couche limite en ce point. 

Solution 



Determinons d’abord la nature du regime d'ecoulement. Pour cela, calculons la 
valeur du nombre de Reynolds au point considere e'est-a-dire a x = 50 cm: 



Re, = 



Ux 0,30x0,50 



= 150000 500000 



v 10 _ " 

L'ecoulement est done Iaminaire et I’epaisseur en ce point s'exprime par I'equation 



3.13: 



S 



5x 

JrT, 



5 x 0,50 
Vl 50000" 



= 0,0064 m 



>l Ces valeurs sont appro ximatives. La valeur transitoire du nombre de Reynolds depend 
de plusieurs facteurs dont le plus significatif est la rugosite du materiau qui favorise un 
passage plus rapide du regime iaminaire au regime turbulent. 
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Exemple 3.5 

De I'eau s'ecoule dans une conduite de section circulaire de diamStre interieur &gal 
A 10 cm. II est demande de determiner la viiesse de I'ecoulement £ la transition. On 
prendra p etm -Id^NsInr et p CM =1000*g7m\ 

Solution 

Pour une conduite de section circulaire, le nombre de Reynolds critique a pour 
valeur2300. Done: 

olJD , Re u. 

R* = - — =* </ = — 

p pD 



d’ou: 



2300 x 10 ~ 
1000 x 0,1 



0,023m/ s 



3.2.5 Nombres sans dimensions 

Les nombres adimensionnels souvent utilises en presence du phenomene de 
convection qu'elle soit naturelle ou forcee sont resumes dans le tableau ci-dessous. 



Nom 


Symbole 


Expression 


Signification des parametres 


Reynolds 


Re 


Forces d' i tier tie UL 

Re - , 

Forces de viscosite v 

Caracterise la nature du regime 
de I'ecoulement (laminaire ou 
turbulent) en convection forcee. 


U: Vitesse caracteristique, 

L : Longueur caracteristique. 

. _ 


Nusselt 




Q echangee par convection 
O echangee par conduction 

hS AT hL 




Prandtl 


— 

Pr 


^ viscosite cinematique 
diffusivite thermhfue 

a k 

Caracterise la distribution des 
vitesses reiativement a celle des 
temperatures e'est-a-dire le milieu 
ou se realise ie transfert. 


a: Coefficient de diffusivite 
thermique. 
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Norn 


Symbole 


Expression 


Signification des paramdtres 


G rash off 



Gr 


^ figP 2 ITAT fig If AT 
Gr - 

P* K- 


/? : Coefficient de dilatation 

1" 1 1 
volumique:^/? - 

g: Acceleration de la 
pesanteur 

AT\ Difference de 
temperature caracteristique. 


Rayleigh 


Ra 


Ra = Pr . Gr 

Remplace le Reynolds en 
convection naturelle. 

— ~ ~ — rr~ 1" 1 


— — rr 1 ^ — ~ 1 — 



3.3 Donnies empiriques de la convection 
3.3.1 Utilisation du nombre de Niisselt 

En etudiant ia couche limite et ses hypotheses, nous avons assume ('existence 
d'une couche mince pres de la paroi dont la vitesse est presque nulle. En s'eloignant de la 
paroi, la vitesse augmente jusqu'a atteindre une valeur constante generalement notee U w . 
En appliquant ces observations qualitatives au transfert de chaleur se produisant entre une 
paroi soiide et un fluide en mouvement, on pourra se faire une idee du profit de la 
temperature. Pres de la paroi, la chaleur ne peut §tre transmise que par conduction; la 
couche au voisinage de celle-ci etant stationnaire et on s'attend a un gradient de 
temperature assez large. En s'eloignant de la paroi, le mouvement du fluide aidera au 
transport de I'energie et le gradient de temperature sera moins important. Pour une couche 
limite turbulente qui s’est form6e pres d'une paroi chauffee, la distribution de ia temperature 
aura la forme representee en figure 3.4. 




paroi 



Figure 3.4: Profil de la temperature 6 I'interieur de la couche limits 
se dSveloppant pr6s d'une paroi soiide 
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Cette discussion suggere une methode devaluation de la quantite de chaleur 
echangee entre la paroi et le fluide. A I'interface (y = 0), le flux de chaleur se transmet par 
conduction et s’exprime done grace a la loi de Fourier (2.5): 



^durtace— jhtide tUridc ^ 



cT 

w,,. 



(3-16) 



Afin de se rapprocher du sens physique du transfert se realisant entre le fluide et la 
paroi, celui-ci est decrit par la loi de Newton (2.24) exprimant un echange par convection. II 
vient: 




(3.17) 



En considerant la temperature a la paroi T p comme temperature de reference, nous 
pouvons ecrire: 



En introduisant une longueur L intervenant comme une caracteristique geometrique 
du corps considere emettant la chaleur, on aura: 



k.tim,i e LS 



( 5T \ 



\ 8 yj v - o 



= hLS(T r -T,): 



hi 






' dT_ 
< £ y)y* o 



T p -T. 

L 






Le rapport — est appele le nombre de Nusselt. Pratiquement, ce nombre permet la 

k f 

determination du coefficient de transfert par convection. En effet. si le nombre de Nusselt 
est connu; h peut etre calcule grace a la simple relation: 



h = Nu-j- 



(3,18) 



3.3.2 Determination du Nusselt 

Sur la figure (3.4), on remarque que le gradient de a temperature est restreint a une 
region mince situee pres de la paroi dont 1‘epaisseur est S T . Afin de simplifier le processus, 
la courbe representant le gradient de la temperature est remplacee par une droite 
(repr^sentee en pointille sur la figure). Cette droite qu est tangente a la courbe represente 
la distribution de la temperature dans une couche fictive depaisseur 6 P qui si elle etait 
stagnante presenterait la meme resistance thermique au flux de chaleur. Dans cette 
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couche, le flux n'est transmis que par conduction. Par unite de surface, celui-ci est exprime 
par la relation: 



O 

s 



= k 



7 -T 

p ** 



fluids 



5 ; 





(3.19) 



Done: 




(3.20) 



Le nombre de Nusselt s'ecrira done: 



Nu 




L_ 

5’r 



(3.21) 



3.4 Methodes devaluation de h 

II existe principalement quatre methodes qui permettent revaluation du coefficient de 
transfert de chaleur par convection h. Toutes les quatre ont contribue a la comprehension 
du phenomene de transfert de chaleur par convection. 

1- une analyse dimensionnelle combinee a I' experl mentation, 

2- une solution exacte des equations du mouvement, 

3- une solution approchee de la couche limite, 

4- une utilisation de I'analogie entre les transferts de chaleur, de 
masse et de moment. 

L'analyse dimensionnelle est mathematiquement simple et a trouve depuis sa 
creation un large champs d'application. Son inconvenient majeur est represente par le fait 
que les resultats qu'elle permet d'obtenir sont incomplets et non utilisables sans I'apport de 
I'experimentation. 

L'analyse mathematique exacte requiert la solution simultanee des equations de 
conservation. La methode presuppose une connaissance precise des mecanismes 
physiques suffisante pour les decrire mathematiquement. La complexity des equations 
rend en fait leur solution difficile. Elle est souvent obtenue grace au concours des 
methodes numeriques. 

Les solutions approchees de la couche limite utilisent une equation unique 
representee par I'integrale des equations du mouvement. Cette methode est relativement 
simple et aboutit a des resultats satisfaisants compatibles avec ceux obtenus grace aux 
methodes exactes. 

L'analogie entre les transferts de chaleur, de masse et de moment est un outil 
generalement utilise pour une analyse en regime turbulent. 



3.5 Analyse dimensionnelle 

L’analyse dimensionnelle est differente des autres methodes par le fait qu'elle 
n’introduit pas d'equations mathematiques a resoudre. Elle permet la combinaison d'un 
certain nombre de variables (ou groupes adimensionnels) qui debouchera sur des relations 
empiriques decrivant des resultats experimentaux d'une maniere acceptable et largement 
utilisable. 
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3.5.1 Determination du nombre de groupes adimensionnels capables de d6crire 
un phenomene: Theoreme de Buckingham ou th&or&me du n 

Afin de determiner combien de nombres adimensionnels sont necessaires afin de 
decrire un phenomene physique, le theoreme de Buckingham peut etre utilise. II definit ce 
nombre (souvent note n d'ou le nom donne au theoreme) comme etant la difference entre 
le nombre total de grandeurs physiques ou parametres intervenant dans la description du 
phenomene {souvent note n) et le nombre des differentes dimensions de base utilisees 
(souvent note m). D’ou: 

7 r =n~tn 



Une valeur nulle ou negative de la difference ( n-m ) signifie que I'ensemble des 
grandeurs physiques envisage ne peut etre ramene a une forme sans dimensions. Pour 
= I'equation de similitude comporte uniquement un seul critere qui sera non 
determinant car ne pouvant etre qu'une constante. Si n-m- 2. I'equation de similitude 
comportera deux nombres sans dimensions (qu'on notera n, et x z ) et pourra s'exprimer: 

f[n x p r,) = 0 ou tt, =i 7 (^- ; ) 

3.5.2 Exemple d'application du theoreme de Buckingham 

A titre d'exemple, essayons de determiner la forme de I'equation de similitude 
permettant la determination du coefficient de transfert de chaleur par convection forcee a 
I'interieur d'une conduite de section circulate. Une certaine experience du phenomene 
considere et de I'utilisation du theoreme du k permet d'assumer que le coefficient de 
transmission h depend principalement de six grandeurs physiques c'est-a-dire que: 

7r(A*,^AA,C*) = 0 (3.22) 

En retenant comme dimensions principals M, L, T et / representant 
respectivement la masse, la longueur, la temperature et le temps, il vient: 

• Diametre du tube: [£>] - L 

• Conductivite thermique du fluide: [£] = ML! t'T 

• Vitesse du flux: [U] - Lit 

• Densite du fluide: [p\ = M / L' 

• Viscosite du fluide: [//] = M / Lt 

• Chaleur massique: [C] = I? ! f 2 T 

• Coefficient de transfert: [h] = Ml t'T 



L'application du theoreme de Buckingham montre un besoin de trois nombres sans 
dimensions pour la description du probleme: 

7 (D,k,U>P,M,C,h)-4 (M,L,T,t) = 3 



Pour Panalyse des dimensions, L’equation (3.22) peut etre ecnte comme: 

n = D°.k b .U s . p d . .h & 



(3.23) 
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En identifiant les dimensions, on aura: 




(3.24) 



Pour que le produit k soit sans dimensions, il est necessaire que ia somme des 
exposants des differentes dimensions soit nulle. soit: 

Pour L: b + d+e + g = 0 
Pour M: a + b + c~3d ~e + 2f = 0 

i ' (3.25) 

Pourt: -3b-c-e-2f-3g=0 V ’ 

^Pour T : -b-f -g = 0 

Ce systeme comporte 4 equations pour six inconnues, II est done possible de 
choisir 3 inconnues pour chaque nombre sans dimensions. La seule restriction dans ce 
choix est que les trois inconnues soient independantes. 



ler choix 

Comme le parametre recherche est le coefficient de convection et afin qu'il 
apparaisse dans le nombre sans dimension, il est choisi g = I. Dans un but de 

simplification des equations, il est pris c = d = 0 . Dans ce cas, la resolution des equations 
(3.25) aboutira a a = 1 , b ~- 1 et e = / = 0 permettant de d^finir le premier nombre sans 
dimensions: 

Nu = A , .Re lf .Pr / (3.26) 

que nous reconnaissons comme etant le nombre de Nusseit. 



26me choix 

h ayant eie defini, il ne devrait plus apparaTtre d’ou le choix de g = 0. En 
choisissant a ~ 1 et / = 0, et en resoivant pour ces valeurs les equations (3.25), il vient 
b- 0, c~d~ 1 et e--\ aboutissant a I’expression du second nombre sans dimensions: 



U Dp 

= 



(3.27) 



que nous reconnaissons comme etant le nombre de Reynolds. 



36me choix 

En prenant cette fois a -c-d - g- 0, b - -1 et e = / = 1 , e'est le 
nombre de Prandtl qui apparaft et qui s’exprime: 

Cp 



(3.28) 
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Par consequent, ['analyse dimensionneile appliquee au cas de la transmission de la 
chaleur lors du mouvement force d'un fluide dans une conduite circulate a permis de 
combiner les six grandeurs physiques dont depend le coefficient de transfert convectif pour 
aboutirfinalementa I'equation de similitude qui s'exprime: 

Nu - K, Re d , Pr 7 (3,29) 

La constante K et les puissances inconnues d et / sont detemninees grace au 
concours de I'experimentation. 



3.6 Formules empiriques couramment utilisees 

Un grand nombre de formules empiriques est disponibie pour determiner le 
coefficient de transmission par convection a travers I'expression du nombre de Nusselt. 
Ces relations dependent notamment du type de convection (forcee ou naturelle) et de la 
nature du regime d'ecoulement du flux (laminaire ou turbulent). L'objet du present 
paragraphe est d’essayer de resumer ces relations et de les exprimer pour le maximum de 
cas possible. 



3.6.1 Convection forcee 

3.6.1. 1 Echange de chaleur le long d'une plaque plane 

Regime laminaire: Re<3.10 ? 

Nu l =0,66 (Re i ) 1 "(Pr) 1 ^ (3.30) 

Regime turbulent: Re^3,10 5 

Ntt, =0,036 (Re j 4 '(Pr) I/3 (3.31) 



Exemple 3.6 

Une plaque mince d'une longueur de 2m et d'une largeur de 1m est sous I'effet d'un 
ecoulement d'air a la vitesse de 1,5m/s et de temperature 20°C dans la direction 
longitudinale. La temperature des surfaces de la plaque est de 90°C. II est demande de 
calculer (a) le coefficient de transmission de la chaleur par convection suivant la longueur 
ainsi que (b) le flux de chaleur transmis par la plaque a I'air. 



Solution 

a/ A 20 ‘C, les caracteristiques de I'air sont: 

p = 1,175 kg i , p = 1,8.1 0”’^/ ms , k = 0,0261^ / mK , C v - 1006// kgK 
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Determinons la nature de i'ecoulement: 

pUL 1,175x1,5x2 



Re 



= 1,96,10 -< 3.10 



-5 



p 1,80,10 - 

II est done laminaire et I'equation (3.30) peut etre appliquee: 

i u n. , , hL k t i is ( pC 

Nu l = 0,66 (Re J (Pr) = — h = 0,66 - (Re J — H 

k L v k 



\l/3 



d'ou: 



0,0260 / j \ i /2 

h = 0,66 - L - — (l, 96.10 ) 



^ 1, 80.1 O^x 1006V ' 3 

v 0,0260 j 



3,36 WhrfK 



b / La chaleur transmise par la plaque (possedant deux parois, Tune superieure et I'autre 
inferieure) a i'air est: 

O = l.h.SAT = 2 x3,36x2,l,5 (90 - 20) = 1411 Watts 



3.6.1. 2 Ecoulement a IMnterieur de tubes cylindriques lisses 



Regime laminaire: Re < 2000 






• D : Diametre interieur du tube, 

• p m , p p - Viscosites dynamiques definies a T m ^T p , 

T +T- 

• T m = P 3 \ Temperature moyenne, 

• Tp. Temperature de la paroi interne du tube, 



0,0668. Re. Pr (DfL) 

Haussen: Nu n =3,66 + ~T 

1 + 0,04 [ Re. Pr. (£»/!)] 




(3.32) 



Sieder et Tate: 



Nv d = 1,86 (Re.Pr) 1 3 




(3.33) 



pour : [Re.Pr.(X>/ Z,)] >-10 

Toutes les proprietes sont definies a T m sauf p p 



Kays: 



Nu d =3,66 + 



0,104.Re.Pr.(Z)/Z) 

1 + 0,016 [Re.Pr.OD/I)] 0 ’ 8 



(3.34) 



pour : [Re.Pr.(Z)/ Z)] -< 100 
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Regime turbulent: Re >- 2000 

Colburn: Nu D = 0,023 (Re) 0 8 (Pr) 1 3 

pour : - L/ D> 60 

- 0,7 < Pr < 100 

- 10 4 ^Re D -< 1,2.1 0 5 

Sieder et Tate : Nv D - 0,023 (Re D ) (Pry-l-^-l 



(3.35) 



(3.36) 



Mc-Adams: 



Nu d = 0,O23.(Re D )' \(Pr)' \ 



f V-14 



[]-(/)/ 1) 0 7 ] 



pour : le regime d'entree dans les tubes. 



(3.37) 



Example 3.7 

Calcuiez le coefficient de transmission de la chaleur par convection ainsi que ie flux 
degage lors de I'ecoulement force d'une huile a la vitesse de 0,5 m/s dans un tube de 
10mm de diametre et de 1m de longueur si ies temperatures moyennes de I'huile et de la 
paroi sont respectivement egales £ 80° C et 20°C. Les caracteristiques de I'huile utilisee & 
la temperature £ laquelle il s'ecoule c'est-a-dire 80 °C sont: p = MAkg hn 1 ; 
ji/ = 30,8.10 ^ kg l ms ; k = 0,1 08 IF/ m°C et C. p - 1 ,%46kJ / kg°C . A la temperature de la 
parol c'estS-dire 20° C, la viscosite de I'huile est: u p = 198,2,10 ~ A kg / ms. 



Solution 

D6temriinons d'abord le regime d'ecoulement: 



Re D 



pUD 

A 



844x0,5x0,01 
30,8.1 0” 4 



= 1370 -< 2300 



Le regime d’ecoulement done laminaire. 



La relation de Haussen (3.32) peut done etre appiiquee. Commengons par caiculer le 
produit: 



Re.Pr.(Z7 / L) = Re 



M C F D 
k L 



1370x30,8.10^x1846x0,01 

0,108.1 



= 721 



0,108 

0,01 



3,66 + 



0,0668.721 
l + 0,04.(72lf 4 



f 30,8 Y ,,iJ 
U98,2> 



/?= ^Nv u 



= 134,6 Wlm 2a C 
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Calcul du flux de chaleur degage: 

O - h.S.&T = h.n.D.LAT = 134,6x^x0,01x1(80-20) = 253,7 Watts 



Remarque 

Le produit [Re. Pr. (£>//.,)] etant superieur a 10, la relation de Sieder et Tate 
representee par I’equation (3.33) peut etre appliquee. Le resultat qu'on trouverait serait 
h = 138 ,%W !m 2o C . L'erreur entre ce resultat et celui trouve grace a I'application de la 
relation de Haussen s'eleve a 3%. 



Exemole 3.8 

On demande de calculer la quantity de chaleur transmise par une eau se d4plagant 
d'une maniere force e dans serpentin constitue d’un tube de 18mm de dia metre. Le debit de 
I'eau est de 0,24kg/s et sa temperature 120°C. La temperature de la paroi interne de la 
conduite dont la longueur est de 3m est consideree constant© et egale a 110 9 C. 

Solution 

Les caracteristiques de I'eau a la temperature de 120°C sont trouvees egales a (cf. 
Annexes): 

p - 945,3 kg! , pi = 2,34.10 ~ 4 kg / ms, k - 0,68 5fV / mK, C p = 4250J / kgK 



Calculons la vitesse de I'eau s'ecoulant dans la conduite et le Reynolds: 
m 0,24x4 

U~ = : 7 = lml s 

pS 945,3 x,;r (0,0 18)' 



Re 



p.U.D 

f-i 



945,3x1x0,018 
2,34.1 0~ 4 



= 7,3.1 0 4 >- 2300 



L'ecoulement est done turbulent. Appliquons la relation de Colburn (equation 3.5). 
Pour cela, il faudra verifier les conditions: 



D 0,018 



= 167 ^ 60 



Virifit 



pi.C p 2,34.1 0“ 4 .4250 

• Pr = = 1,45 ^ 100 et >- 0,7 

• Veripe 



k 0,685 

Re = 7,3.1 0 4 ^ 1,2,1 0 5 et ^10 4 



Verifit 



Done: 

h=~ > Nu D = "q^Y^O, 023.(7, 3.10 4 ) ia .(l,45) 1 3 = 7702W lm 2 K 
O = h.n.D.LAT = 7702.^.0,018.3.(120-110) = 13000^ = 13 kW 
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Remarque 

La relation de Sieder etTate (equation 3.36) permet aussi de calculer le nombre de 
Nusselt et par consequent le coefficient de transmission et le fiux. Son application aboutit 
aux resultats suivants (jj = 2,57.10'* kg /ms): 

Nu d = 212,5 , h = 8086,8 ^ IntK , <J>=\\lkW 

L'erreur entre les resultats fournis par ies deux relations (Colburn et Sieder & Tate) 
est de 5% dans ce cas. 



3. 6.1. 3 Ecoulement dans les espaces annulaires 




Figure 3.7: Ecoulement dans I'espace anruiaire forme par les deux conduites 



Nu dh = 0,023 (Re) (Pr)" 



avec: • Re DH = 



II ..DH 



et 



Nu dh - 



h.DH 



v "" k 

• DH: Dia metre hydra ulique. Dans ce cas: DH = D 2 - D s 

• n - 0,4 pour chauffage (j. >- rj 

• n = 0,3 pour refroidissement 



(3.38) 



3.6.1 .4 Ecoulement perpendiculaire & un tube 

Hilpert : Nu D = C.(Re D ) 



(3.39) 



R e D 


C 


m 


1-4 


0,891 


0,330 


4-40 


0,821 


0,385 


40-4000 


0,615 


0,466 


4000-40000 


0,174 


0,618 


40000-250000 


0,0239 


0,805 




Figure 3,8: Ecoulement perpendiculaire 
A une conduite 
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3.6.1. 5 Ecoulement perpendiculaire a une rang£e de tubes 

oloTo o oTo 

|OjO 0|0 op 

Figure 3.9: Ecoulement perpendiculaire £ une rangee de conduites 



Colburn: Nv D = 0,33 (Re p f’ 6 (Pr) 1 



(3.40) 



Exemole 3.9 

Une barre de sect/on circulaire de diametre 15mm est refroidie par un couranf d’air 
transversal de vitesse Im/s et de temperature 20° C. Calcuiez la quantity de chaleur 
transmise a I'air par unite de longueur de la barre si la temperature de la paroi de cette 
dernidre est de 80°C. 




Solution 

A 20‘ C, les caracteristiques de I'air sont (cf. annexes): 

p ~ 1,2 kg l m 3 , p - 1,8.10"' kg I ms, k - 0,0259FF / mK, C p = 1006J / kgK 



L'ecoulement etant perpendiculaire a la barre, la relation d'Hilpert representee par 
I’equation (3.39) peut etre appliquee. Calculons d'abord le nombre de Reynolds: 



Re D 

d'ou: 



P-U.D _ l 
M 



h = ~Nu D 



20x1x0,015 
1,8.10 
k 



,-5 



D 



c(Re„) 



= 1000 => 

w 0,0259 
0,015 



C = 0,615 et/n = 0,466 
0,615(1000) ( ' ,4ti6 =26,55 WfrtrK 



Finalement: 

<D =A.5.A7'= 26,55xn-.0,015(80-20) = 75 W/m 



3.6.2 Convection naturelle ou fibre sur parois isothermes 

En paragraphe 3.5.2, il a ete demontre qu'une correlation de trois nombres 
adimensionnels (en (’occurrence: Re, Pr et Nil) est necessaire afin de decrire un 
phenomene de transfert de chaleur par convection forcee. Dans le cas de la convection 
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naturelle ou libre, il peut etre demontre qu'une correlation entre nombres adimensionnels 
est aussi indispensable. Dans ce cas, les nombres adimensionnels sont ceux de Rayleigh, 
Grashoff et Prandtl. On notera que le nombre de Rayleigh a remplace celui de Reynolds 
pour la determination du regime d'ecoulement (laminaire ou turbulent, cf. Tableau de la 
page suivante), 

Le present paragraphs propose les relations pouvant decrire un phenomene de 
transfert de chaleur par convection naturelle sur des parois isothermes c’est-a-dire dont les 

surfaces sont a une temperature constante (T p = Cs1e). 



Pour un aaz auelconaue; 


Nu f 


= A.(Gr f .Pu f Y 


(3.41) 


Dans le cas de I'air: 


Ntfj 


= B- (Cr,)' 


(3.42) 



Re marques 

1- Le nombre de Prandtl pour Pair est generalement considers constant et sa valeur est 
prise egale a 0.7. 

2- L'indice '/’ indique que Ea temperature du fluide est prise egale a celle du film situe 

T + T„, 

pres de la paroi c'est-a-dire que T f = — . 



Geometrie 


Ra f 


B 






Obs. 


Plans et Cylindres 
verticaux 

Plans: L=Hauteur 
Cvlindres: L=Lonoueur 


Laminaire: 1 0 4 — 1 0 9 
Turbulent: 10 9 -10 lJ 


0,59 

0,13 






Gr evalue 
sur L 


Cylindres 
horizontaux 
de diam6tre: D 


Laminaire: 10‘ ! - 10 9 
Turbulent: 10 9 - 10 u 








Gr evalue 
sur D 


Surfaces planes 
horizontales de 
longueur L dans le 
sens de l'6coulement 


Laminaire: 10 5 -2.10 7 
Turbulent: 2.1 0 7 - 3.1 0 ,fl 

- 


0,59 

0,13 


0,25 

0,33 


0,54 

0,12 


Gr lvalue 
sur L 
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Exemole 3.10 

Calculez ies pertes de chaleur par unite de temps et par metre carte de surface d'un 
6changeur horizontal dont le corps cylindhque est refroidi par un courant d'air Iibre de 
temperature 30“C. Le diam&tre exteme de la conduite est de 400mm et la temperature de 
sa paroi 200°C, 



Solution 



On est dans le cas d’une convection Iibre avec comme fluide de I'air. La relation 
(3.42) peut done etre utilisee . Le choix des constantes B et m necessite la connaissance 
du regime d’^coulement qui dans le cas d’un d’un echange par convection naturelle est 
decrit par le nombre de Rayleigh. Caiculons d'abord la temperature du film de fluide: 




T, + T 



paroi 



2 



200 + 30 
2 



= 115°C= 388AT 



A cette temperature, Ies caracteristiques de I'air sont (cf. Annexes): 



p - 0,885&g/ m i , p- 2,2.1 O' '' kg f ms, k = 0,0331^ / mK, C p - 1 013-7 / kgK 



Done; Gr, - 



J - 88 ^1.(0,885)+(0,4y .(200-115) 



( 2 , 2 . 10-7 



= 2,22.10 



Pr , 



M-C p 2,2.1 O' 5 .1013 
k ~ 0,033 



= 0,675 



d’ou: 

Ra f = Gr f . Pr,. = 2,22.10 s .0,675 = 1,5.10 s + 10 9 



Le regime est done laminaire. La conduite etant horizontale, les constantes de 
I’equation (3.42) sont done egales a (cf. Tableau du paragraphe 3.6.2): 

B = 0,49 et m = 0,25 



Finalement: 



h 





0,033 

0,4 



.0,49.(2,22.1 0 s ) 



(i_2 5 



4,93 W!m 2 K 



Et: <&=hAT= 4,93.(200-30) = 838 W / m 2 



Exemole 3.11 

Une plaque horizontale de (2x3)m2 de surface est 0 hentee vers le haut par sa paroi 
emettrice de chaleur. Calculez le coefficient de transmission de la chaleur a I'air ambient 
calme sachant que la temperature de la plaque est 120°C et celle de l‘air34°C. 
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3.6: 

La temperature en tout point d’une mince plaque plane placee paralieiement a un 
courant d'air est egale £ 365K. La Vitesse du courant d'air principal est de 60m/s et sa 
temperature 270K. La plaque possede une langeur de 0,6m et une longueur de 0,45m 
dans la direction de I'ecoulement. En nbgligeant les effets de bond de la plaque et en 
assumant que I'ecoulement dans la couche llmite change brusquement du laminaire au 
turbulent 3 un Reynolds de 4. 10' , trouvez: 

1- la distance £ laquelle se fait la transition du laminaire au turbulent , 

2- les coefficients de transfert convectif moyens (en tenant compte du nombre de Reynolds 
moyen) dans les regions laminaire et turbulente, 

3- le flux de chaieur de la plaque entiire considbrant les deux faces. 

Donnbes: On prendra les caracteristiques de I’air comme: 

[p = 1,2kg I w 1 , C-p = 1010.// kgK , Pr = 0,71 
l* = 0,024 W/mK , // = 1,74.10 '’kg! ms 



R6p.: 1- L 

2- h r 



transition 



= 0,097m 



= 65,15 W ! nr K , h T 



3- <D = 6,8 kW 



l$\,9W!m l K 



3.7: 

Un radiateur est assimile a une plaque verticale de (Imxlm). Que devrait etre sa 
temperature pour dissiper 0,5kW dans un ambiant d'air de temperature 20°C s'ecoulant par 
convection forc&e £ la vitesse de Im/s ? 

R6p.: T - 149’C 



3.8: 

Les deux equations suivantes ont ete propos&es par Hansen en 1943. La premiere 
est applicable 6 des cas d'echange de chaieur par convection fovc&e dans un 6coulement 
laminaire d i’interieur de conduites cylindriques, tandis que la seconde s'applique pour des 
nombres de Reynolds se situant dans la region de transition c'est-a-dire 2000 -< Re ■< 8000 
ainsi que pour des Reynolds plus importants. 

Comparez les valeurs du nombre de Nusselt predits par les deux Equations pour 

- Re=1000 Pr=1 pour (D/L)=0,08 et (D/L)=2 

- Re=3000 Pr=1 pour (D/L)=0,1 

~Re=20000 Pr=1 pour (D/L)=0,01 

avec celles obtenues paries relations empiriques appropriees. 



(V 

( 2 ) 



Nu = 3,65 + 



0,668(Z)//,)RePr 
l + 0,04[(ZJ/I)RePrf /3 



Nu = 0,1 06 [(Re) 2 3 - 125](Pr) m [l + (Di L ) 2/3 ] [J^ 
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On assumera que le fluide est de I'eau s'ecoulant a la temperature de 288K a 
I'interieur d'une conduite de temperature 336K. On assumera aussi que sa viscosite 
dynamique vane avec la temperature suivant le tableau suivant: 



T(°C) 


p (kg/ms) 


0 


1,78 . 10~* 


10 


1,00.10 


20 


1,00 . 10 3 


40 


0,651 . 10~~ 


60 


0,469 . 10- i 


80 


0,354 . 10 - 3 


100 


0,281 . 10~* 



Rep.: 



Re = 1000 

• 1 => Re.Pr.(Z) / L) = 80 -< 100 



[D/L = 0,08 



Re = 1 000 

• I „ => Re.Pr.(£>/£) = 2000^ 10 

\D!L = 2 





= 9,1 




NUswder&Tau 


8,8 


Nu (]) = 


7 




N U S)e&r&Tate ~ 


25,78 




23,5 





f Re = 3000 
[D/I = 0,1 



=> Re -< 1 0000 



{ ^ f Sieder±TMe ~ 

Nu (2) = 12.8 



15.30 



Re = 20000 

< , => Re > 1 0000 etL / D >~ 60 

[D/L = 0,01 



f Nucm™ = 6 ^ 4 
[Nu^ = 81,7 



3 . 9 : 

Calculez par deux formules differentes la quantite de chaleur transmise par I'eau 
s'ecoulant d'une maniere forcee dans un serpentin constitue d'un tube de 18mm de 
diametre. Le debit de I'eau est de 0,24kg/s et sa temperature 120°C. La temperature de la 
paroi interne de la conduite dont la longueur est de 3m est considerPe constante et egale e 
1 10°C. Comparez les resultats obtenus et expliquez la difference qui les separe. 

Les caracteristiques de I'eau sont: 



A 100°C 



A 120‘C: 



[C, = 4216/ kgK ; k = 0,6801*' / mK 
= 960,6£g / m' ; p- 0,28 1.1 0 kg / ms 

| C p = 4250/ kgK , k - 0,6851*' / mK 
Ip = 945,3kg! m i ; p = 0,234.10 *kg / ms 



R&LI *^ = 12,9 kW ; ® Siedtr&TM = UJkW 

L’erreur de 6,2% trouvee est principalement due a I'empirisme des 
relations utilisees. 
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3 . 10 : 

De fair d une pression de 2 atmospheres et a une temperature de 200 °C s'6chauffe 
lorsqu'il passe dans une conduite de 1,58cm de diam&re a une vitesse de 15m/s. 

1- Calculez la quantity de chaleur transmise par unite de longueur de la conduite si le flux 
thermique est maintenu constant & la paroi dont la temperature est de 20°C au-dessus de 
celle de fair tout le long de la conduite, 

2- Quelle serait I'augmentation de la temperature sur une longueur de 3m si on garde les 
memes conditions ? 

Indications: 1 - L 'air est a ssimite £ un gaz parfait (P = p.R.T) 

2- L'energie peut s’exphmer pan <b = m.C p ,AT 

R6p.: 1- O=l00Wfm ; 2-AT=67°C 



3 , 11 ; 

1- Un tube horizontal de 70mm de diametre traverse une enceinte ferm&e de 5m de 
longueur remplie d’air d la temperature de 20°C. La temperature de la parol de la conduite 
est de 100 a C. Calculez le flux 6chang6 par la conduite par convection. 

2- En ouvrant les fenetres de I’enceinte, un courant d’air d’une vitesse de 2m/s 
prend place. En assumant que les temperatures de fair et de la conduite varient en 
consequence pour devenir eg ales a 10°C et 90°C respectivement, calculez dans ce cas le 
flux echang6. 





10°C 


20°C 


50° C 


60°C 


90°C 


(*„? / W 3 ) 

Mnr 

QV (Jikg.K) 

k mr {Wim.K) 


1,24 

1,76.10~- 

1006 

0,024 


1,20 

1,81.10 “ 5 
1007 
0,026 


1,09 

1,94.10 * 
1008 
0,028 


1,03 

2,02.10' 

1009 

0,030 


0,99 

2,08. 10~ % 
1010 
0,032 



Rep.: 1-$ = 703,31fF ; 2- O = 876 W 



3 . 12 : 




I 



e=0,01m 



1- Soitle murreptesente ci-dessus. En assumant une conduction unidimensionnelle 
et co nnaissa nt les valeurs des temperatures aux parois T pl et T ?l ainsi que les 
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conductibilitps thermiques de ses differentes sections, calculez le flux de chaleur transmis d 
t ravers le mur. 

2- En assumant I'Gpaisseur du mur negligeable devant ses autres dimensions , celui- 
ci est assimile a une plaque horizontale au-dessus de iaquelle s'ecoule librement de I’air & 
la temperature de 1000K dans le sens de la longueur. La temperature de la parol de la 
plaque est de 1074K. On demande de determiner le flux de chaleur echange entre I'air et 
la plaque . 

3- Trouvant que le flux echange par convection fibre est insuffisant, un ventilateur est place 
au-dessus de la plaque permettant de faire acquerir £ I'air une vltesse de 20m/s. Les 
temperatures de I'air environnant et de la plaque diminuent en consequence et deviennent 
6gales a 890K et 980K respectivement. On demande de caicuier dans ce cas le nouveau 
flux de chaleur 6changP entre I'air et la plaque. 

Donnees: 

T pl = 1850°C ; J;,=65°C 

\k A - 180 W ! mK ; k B =4SWimK ; k c =90WfmK 
k D =60 WhnK ; k E =\20 W!mK ; k F =\65WlmK 

R6p.: 1- 1400 W 2-59W 3- 5760JF 



3.13 : 

1- La temperature d'une plaque verticale parallelement a Iaquelle s'Pcoule un gaz dont les 
caracteristiques sont resumSes dans le tableau est de 385K, la temperature du gaz etant 
31 5K. La plaque possede une largeur de 0,1m et une hauteur de 0,2m dans la direction 
de I'ecoulement. II est demande de determiner 

la- la nature de I'ecoulement, 

1b- le coefficient de transmission par convection, 

1c- le flux de chaleur dSgage par fa plaque. 

2- On place un ventilateur au-dessus de la plaque de telle fagon que le gaz environnant 
acquiert une vitesse de 60m/s. En assumant que les temperatures de la plaque et du 
fluide diminuent en consequence et deviennent egales respectivement d 340K et 260K, 
dSterminez: 

2a- la nature de I'ecoulement, 

2b- le coefficient de transmission par convection, 

2c- le flux de chaleur degage par la plaque. 

3- Quelle conclusion pouvez-vous tirerdes resultats obtenus ? 

Donnees: caracteristiques du gaz: 



T(K) 


p (kg ! rrri ) 


M (kg / m. .v ) 


C, (Jlkg.K) 


k{Wim.K) 


250 


2,166 


12,6 10' s 


803,9 


0,0129 


300 


1,797 


15,0 10 _e 


870,9 


0,0166 


350 


1,536 


17,2 10" 6 


900,2 


0,0205 


400 


1,342 


19,3 10^ 


942,0 


0,0246 
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Ftep.: la- laminaire ; 1b- h = 5 WlnrK ; 1c-<$=7W 

2a- Turbulent ; 2b- h = 233 WlnrK ; 2c- 372,8 W 

3- La convection forcee augmente d'une maniere sensible le flux 
de chaleur degage par la plaque 



3.14 : 

On sait que les echanges par convection libra sont d&cnts par une relation du type 
Nu=f(Gr,Pr). Prenons ie cas de fair pour Sequel: 

p = 1,25.10 ' J CGS a Q°C ; p = 1, 74.1 O ' 4 CGS a 0°C 

k=0,58 calories CGS a 0°C ; Pr=0,71 



Gr = 



fi.g-AT.p--L 3 

M 2 



L: Dimension lineaire de I'objet considere. 



log l0 (Nu) 


0,16 


0,28 


0,42 


0,59 


0,78 


1,02 


log Kl (C/\Pr) 


0 


1 


2 


3 


4 


5 




1,27 


1,52 


1,77 


2,02 


2,33 


2,70 


6 


7 


8 


9 


10 


11 



1- En prenant L=3m, donnez pour ST =5°C, 10°C et 15°C ia vaieurdu NOsselt Calcuiez ie 
coefficient de transmission h enW i cm 2o C et tracez la courbe h = f (ST) . 

2- Un mur est constituG d'une paroi en verre de 5mm d'epaisseur et de conductivity 
thermique k- 10~ 2 W/ cm 0 C , Calcuiez sa resistance thermique en W i CGS . 

3 - Calcuiez, pour un fiux de chaleur de 20W/mr, 40W / m 2 et 60W ! m 1 la difference de 
temperature entre les deux feces de la plaque de verre et en assumant que la valeur de h 
est la m$me des deux cotes tracer la courbe repnesentant la variation du flux en fonction de 
la difference de temperature. 

4 - Remplacez le simple vitrage par deux vitres identiques separees par un vide d'air de 1cm 
d'epaisseur. Tracez la courbe repnesentant la variation du flux de chaleur par unite de 
surface en fonction de la difference de temperature. 

R6p.: 1- Gr = 5.10 7 .Ar ; 2- R = 50° C / W 

3- ST = 0,1/ 0,2/ 0,3 ; 4- ST = 4266,7$ 



3.15 : 

Une conduite honzontale de 30cm de diametre, de 3m de longueur et de 250“ C de 
temperature se trouve dans une piece ou regne une temperature de 15°C. Calcuiez la 
puissance echangee par convection fibre entre la conduite et fair de la piece. 



Ftep.; $ = 4,45 k-W 
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3 . 16 : 

Le mur d’une maison a 6m de hauteur et 10m de longueur. Sous I'&chauffement du 
soleil, sa temperature atteint 40°C. La temperature ambiante est de 20°C. Calculez le flux 
^change par convection naturelle par le mur 

R6p.: cd = 4,1 kW 



3 . 17 : 

Calculez le coefficient de transmission de la chaleur d’une plaque verticale de 
hauteur 2m dont la temperature de la surface est de 100°C si elle 6tait laiss&e A fair fibre 
ayant pour temperature 20*C. 



